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SECTION  I 
INTRODUCTION 


The  results  for  this  research  program  are  summarized  in  this  Final 
Technical  Report.  The  primary  goal  of  this  research  effort  was  to  develop 
an  efficient,  flexible  computer  algorithm  for  computing  the  coupling 
between  pairs  of  Army  antennas.  Emphasis  was  placed  on  the  development  of 
a  practical  algorithm  for  analyzing  the  coupling  between  pairs  of  directive 
antennas  when  both  antennas  are  operating  within  their  design  frequency 
bands  and  are  located  within  the  near  field  of  each  other.  The  effects  of 
scattering  from  nearby  support  structures  and  other  antennas  on  antenna 
coupling  were  also  addressed  in  this  research  program.  Clearly,  the 
electromagnetic  coupling  analysis  of  Army  antennas  must  be  based  on 
analytical  techniques  that  are  efficient  and  accurate  for  analyzing  near¬ 
field  coupling  and  scattering  phenomena.  The  analytical  techniques  should 
also  be  flexible  enough  to  permit  extension  to  handle  coupling  analysis 
over  wide  frequency  bands  at  both  in-band  and  out-of-band  frequencies. 
Although  major  emphasis  was  given  to  coupling  between  directive  antennas 
unde;  both  "clear-site"  conditions  and  "partially-obstructed"  conditions, 
the  coupling  model  described  herein  is  also  applicable  to  antennas  of 
moderate  or  low  directivity. 

One  can  readily  identify  several .basic  methods  that  have  been  applied 
to  the  analysis  of  antenna  coupling  for  restricted  classes  of  antennas. 

These  are  the  Method  of  Moments  (MOM)  [1,2]  ,  Aperture  Integration  (Al) 

[3,4],  the  Geometrical  Theory  of  Diffraction  (GTD)  [5-7],  and  the  Plane 
Wave  Spectrum  (PWS)  analysis  [8-11].  Methods  based  on  a  spherical  harmonic 
wave  expansions  have  also  been  formulated  for  computing  clear-site  coupling 
[10,11].  Finally,  a  method  based  on  the  Spherical  Angular  Functions  (SAF) 
for  the  antennas  and  scattering  obstacles  was  developed  on  this  contract  [12] 

The  MOM  is  efficient  for  computation  of  coupling  between  electrically- 
small  (low  directivity)  antennas  operating  on  a  clear-site  or  operating  in 
the  presence  of  electrically-small  scattering  obstacles.  The  MOM  cannot  be 
economically  applied  to  antennas  of  moderate  or  high  directivity  due  to  the 
excessive  number  of  wire  elements  required  for  accurate  modeling.  The  AI 


and  GTD  methods  have  been  applied  in  the  analysis  of  electrically-large 
(directive)  antennas,  and  the  GTD  has  also  been  applied  to  electrically- 
small  antennas  mounted  on  vehicles  and  aircraft.  However,  the  AI  and  GTD 
methods  require  recomputation  of  the  electric  fields  in  the  near-field 
whenever  either  antenna  is  rotated  or  when  either  antenna  or  obstacle 
locations  are  changed. 

The  PWS  analysis  overcomes  many  of  the  disadvantages  of  the  AI  and  the 
GTD  analyses  for  directive  antennas  whose  locations/orientations  are 
varied.  In  general,  spectral  analysis  techniques  are  particularly  well 
suited  for  efficient  and  accurate  characterization  of  a  wide  variety  of 
near-field  coupling  and  scattering  phenomena  involving  electrically-large 
antennas  and/or  scattering  objects.  To  date,  the  PWS  analysis  has  been 
employed  to  analyze  the  effects  of  near-field  scattering  obstacles  on 
antenna  performance  for  a  wide  variety  of  situations  that  would  be 
comparatively  awkward  and  less  efficient  to  analyze  via  the  Method  of 
Moments  (MOM),  aperture  integration  and/or  the  Geometrical  Theory  of 
Diffraction  (GTD).  The  PWS  analysis  has  also  been  applied  to  analyze  the 
coupling  between  a  pair  of  near-field  directive  antennas,  under  both 
"dear-site"  conditions  [8-11]  and  with  an  intervening  metal  cylinder 
located  between  them  [10].  The  success  of  the  PWS  analysis  for  numerical 
analysis  of  near-field  coupling  and  scattering  for  directive  antennas  and 
scattering  obstacles  is  due  largely  to  the  speed  of  the  Fast  Fourier 
Transform  (FFT)  that  is  employed  when  the  integrations  over  the  transverse 
wavenumbers  kt  are  confined  to  "visible  space"  limits  | kt|  <_  k0  (k0  * 

2tt/A,  where  A  is  the  wavelength)  [13]. 

It  is  important  to  distinguish  between  the  rigorous  Plane  Wave 
Spectrum  analysis  presented  by  Kerns  [14]  and  the  approximate  analysis  that 
is  used  in  practice.  In  practice,  researchers  typically  work  only  with  the 
PWS  over  a  portion  of  the  "forward"  hemisphere  and  ignore  the  propagating 
electromagnetic  fields  at  wide  angles  in  the  front  hemisphere  and  the 
"evanescent"  fields  in  the  "invisible"  space.  This  procedure  has  proven  to 
give  accurate  results  when  the  coupling  or  scattering  integral  over  the 
excluded  energy  is  small  compared  to  the  corresponding  integrals  over  the 
included  energy,  as  is  the  case  for  many  coupling  or  scattering  situations 
involving  directive  antennas  and  3catterers. 
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There  are  some  potential  difficulties  associated  with  the  approximate 
PWS  analysis  for  certain  commonly-encountered  orientations  of  closely- 
spaced  antennas/scatterers  for  situations  where  one  of  the  antennas,  both 
of  the  antennas,  or  one  or  more  of  the  scattering  obstacles  have  low 
directivity.  Consider  first  the  situation  depicted  in  Figure  1-1  for  two 
directive  antennas.  The  PWS  analysis  is  readily  applied  to  the  coupling 
situation  depicted  in  Figure  1-1  for  both  the  clear  site  situation  shown  in 
the  figure,  as  well  as  with  a  directive  scattering  obstacle  located 
anywhere  between  the  two  antennas,  provided  the  rotation  angles  are  not 
"too  large".  When  the  antennas  are  closely-spaced  and  the  peak  of  the 
mainbeam  approaches  to  within  about  2  or  3,  3-dB  beamwidths  of  90  degrees, 
it  may  be  necessary  to  extend  the  PWS  coupling  integral  integration  limits 
beyond  the  usual  "visible  space"  limits  |kc|-k0.  Although  this  procedure 
is  analytically  straightforward,  it  can  be  inconvenient  numerically  because 
it  would  be  necessary  (1)  to  compute  the  PWS  spectrum  functions  for  |kt|  > 
kQ  and  (2)  to  compute  the  coupling  integral  which  has  a  numerical 
singularity  in  the  integrand  at  jkt|  *  kQ.  Alternatively,  one  could 
possibly  choose  a  new  plane  centered  on  the  line  joining  the  centers  of  the 
two  antennas,  or  perhaps  two  or  more  "canted"  planes  which  allow  one  to 
avoid  using  portions  of  k-space  for  |kt|  >_  kQ.  However,  values  of  coupling 
computed  numerically  for  closely  spaced  antennas  may  be  sensitive  to  the 
choice  of  canted  planes  and  to  numerical  "filters"  that  are  employed, 
deliberately  or  inadvertently,  in  the  computer  processing.  The  practical 
concerns  just  described  are  even  more  cogent  for  arbitrary  coupling 
situations  such  as  the  one  shown  in  Figure  1-2  for  two  directive  antennas, 
and  for  antennas  of  low  directivity  in  general.  In  fact,  the  coupling 
situation  depicted  in  Figure  1-2  probably  occurs  more  frequently  at  Array 
field  sites  than  does  the  one  depicted  in  Figure  1-1,  since  care  is 
ordinarily  taken  to  avoid  pointing  the  mainbeam  of  an  antenna  in  the 
general  direction  toward  a  susceptible  antenna. 

In  light  of  the  foregoing  considerations,  an  alternative  analysis 
referred  to  as  the  Spherical  Angular  Function  (SAF)  analysis  has 
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been  formulated  to  accurately  compute  coupling  data  for  arbitrary  antenna 
arrangements . 

Based  on  the  technical  considerations  mentioned  in  the  preceeding 
paragraphs ,  the  approach  that  was  taken  to  achieve  the  overall  objective 
stated  in  the  first  paragraph  consisted  of  the  following  two  major  tasks: 

TASK  1:  Theoretical  and  numerical  investigations  of  the  Spherical  Angular 

Function  coupling  analysis  were  performed  (1)  to  develop  the  theory  and 
equations  for  analysis  of  arbitrary  antenna  arrangements  and  (2)  to  derive 
suitable  approximate  coupling  equations  for  EMC  engineering  applications 
that  employ  the  use  of  selected  subsets  of  the  full  spherical  angular 
spectrum. 

TASK  2:  Based  on  the  results  of  Tasks  1  and  2,  a  practical,  user-oriented 
antenna  coupling  algorithm  for  application  to  Army  antennas  was  developed 
and  implemented. 

The  theory  and  equations  for  the  rigorous  SAP  are  presented  and 
discussed  in  Section  II.  Useful  approximate  analyses  and  equations  are 
presented  in  Section  III.  Illustrative  antenna  coupling  curves  computed 
via  computer  algorithm  GTSAF  and  auxiliary  programs  are  presented  and 
discussed  in  Section  IV.  Concluding  remarks  and  recommendations  are 
contained  in  Section  V. 
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SECTION  II 

SPHERICAL  ANGULAR  FUNCTION  ANTENNA  ANALYSIS 


A»  Introduction 

A  coupling  analysis  is  sought  which  is  readily  applied  to  arbitrary 
orientations  and  arrangements  of  antennas  located  within  the  near-field  of 
each  other.  No  restrictions  are  made  concerning  the  directivity  of  the 
antennas.  Multiple  scattering  among  the  antennas  is  ignored  in  the 
analysis  presented  herein.  The  coupling  analysis  is  first  conduced  for 
antennas  that  are  operating  on  a  "clear  site"  free  from  scattering 
structures.  The  coupling  analysis  is  then  extended  to  handle  effects  of 
scattering  from  obstacles  that  are  located  within  the  near-field  of  the 
antennas . 

The  coupling  analysis  described  herein  employs  the  Spherical  Angular 
Functions  (SAF)  of  each  antenna.  The  SAF,  denoted  as  F  ( 9,  <P  ) ,  is  related 
to  the  far-field  electric  field  E^f(0,$)  by  the  equation 

iff(ef*  )  -  — —  '<«.♦  >  (2-D 

r 

where  (r ,  0,  <j>  )  are  the  range,  elevation  angle,  and  azimuth  angle  of  a 
selected  point  in  the  far-field  region  of  the  antenna.  Hence,  F(  9,  <f>)  is 
then 

F(  9,  4)  -  re^  If  f  (  9,  «  (2-2) 

where  we  note  that  E ^ ^ (  0 ,  <{>)  and,  hence,  F(  9,  <J>)  are  analytic  over  the  entire 
far-field  sphere  for  any  realizable  antenna.  This  is  a  desirable  property 
since  the  analysis  is  based  on  a  coupling  integral  involving  a  convolution 
of  the  SAF's  for  the  two  antennas  integrated  over  the  entire  far-field 
sphere.  The  theory  and  equations  for  the  SAF  analysis  are  presented  and 
discussed  in  Subsections  B  and  C  for  antennas  operating  on  a  clear  site. 

The  theory  and  equations  for  antennas  operating  in  the  presence  of  near¬ 
field  obstacles  are  sunmarized  in  Subsection  D. 

B,  Coordinate  Space  Coupling  Integral 

A  convenient  starting  point  for  the  analysis  is  the  "reaction",  or 


"coupling”,  integral  due  to  Runsey  [  15  ].  The  coupling  integral,  denoted  as 
((b»a))  ,  is  equal  to  the  open  circuit  voltage  Vba  induced  in  the 
"terminals"  of  the  receiving  Antenna  B  by  the  transmitting  Antenna  A.  The 
terminals  referred  to  here  can  be,  for  example,  the  gap  between  the  inner 
and  outer  conductors  of  a  coax-to-waveguide  feed,  or  a  coupling  iris,  or  a 
selected  "terminal  plane"  inside  the  waveguide.  The  coupling  integral 
given  as  Equation  49  of  Reference  15  is  written  here  as  the  surface 
integral 


Vba  -  (( b,a))  - 

-  J  f[  Ea(R+?')  x  Hb(f ’ )  +  Ha(R+? ' )  x  Eb(r’)]*  n’  ds'] 

*S'  (2-3) 

where  S'  is  an  arbitrary  surface  enclosing  Antenna  B  and  n'  is  a  unit 
vector  normal  to  the  surface  and  pointing  outward,  as  shown  in  Figure  2-1. 

Ea  and  Ha  are  the  electric  field  and  the  magnetic  field  "intensity" 
produced  by  Antenna  A  on  the  surface  S'.  The  interpretation  of  Eb  and  Hb 
depends  on  whether  the  receiving  antenna,  Antenna  B,  contains  gyrotropic 
material  [16]. 

If  Antenna  B  does  not  contain  any  gyrotropic  material,  Eb  and  Hb  are 
the  fields  that  Antenna  B  would  establish  on  the  surface  S'  when  used  as  a 
transmitter  with  unit  input  current.  In  this  case,  the  transmitting  and 
receiving  patterns  are  said  to  be  equal.  More  precisely,  the  transmitted 
complex  fields  Eb  and  Hb  established  at  a  point  are  proportional  to  the 
complex  voltage  response  of  the  antenna  when  used  as  a  receiver  with  an 
idealized  infinitesimal  dipole  source  located  at  If  Antenna  B  contains 

gyrotropic  material,  such  as  a  ferrite  device,  Eb  and  Hb  are  the  fields 
established  on  the  surface  S'  with  a  unit  input  current  and  with  the 
direction  of  the  applied  magnetic  field  reversed  from  the  direction 
normally  used  during  transmission.  In  this  case,  the  receiving  pattern  is 
said  to  be  equal  to  the  "transpose"  transmitting  pattern.  More  precisely, 
the  transmitted  fields  Eb  and  Hb  es"ablished  when  the  biasing  magnetic 
field  is  reversed  through  the  ferrite  are  proportional  to  the  voltage 
response  when  the  antenna  is  used  as  a  receiver  with  idealized 
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Sketch  depicting  antenna 


infinitesimal  dipole  source  located  at  T'  and  the  biasing  nagnetic  field 
applied  in  the  direction  for  normal  operation. 

The  polarization  states  of  the  two  antennas  must  be  specified  with 
respect  to  a  chosen  coordinate  system  in  order  to  perform  the  vector 
operations  indicated  in  Che  integrand  of  the  coupling  integral.  If  Che 
polarization  of  each  antenna  is  resolved  in  cartesian  components  with 
respect  to  non-rotating,  parallel  cartesian  coordinates  fixed  at  the  center 
of  rotation  of  each  antenna,  the  vector  operations  will  be  independent  of 
the  relative  location  of  the  antennas.  However,  the  vector  operations  will 
be  functions  of  relative  rotation  angles  of  the  antennas.  If  the  antenna 
polarizations  are  each  resolved  in  terms  of  spherical  components  with 
respect  to  the  parallel,  non-rotating  cartesian  coordinate  system  at  the 
centers  of  rotation  of  the  two  antennas ,  the  vector  operations  will  depend 
on  both  the  relative  location  and  the  antenna  rotation  angles.  However, 
the  vector  operations  can  be  readily  performed  via  spherical  transformation 
and  rotation  matrices.  The  computer  program  will  be  designed  to  compute 
coupling  versus  relative  orientation  angles  for  a  fixed  relative  location 
of  the  two  antennas.  Consequently,  spherical  polarization  components  are 
selected. 

Antenna  patterns  are  customarily  recorded  and/or  computed  for  (9,tf>) 
polarizations.  The  r-polarlzed  components  of  £  and  ?  can  be  expressed  in 
terms  of  the  (9,$)  components  via  the  Maxwell  curl  relations 


E 

H 


— —  (V  x  H),  and 
ju>e 


1 

jwU 


(7  x  E) 


which  lead  to  the  equations 


E  r 
r 


1 

jwe 


rsin9 


39 


(H  sin(0)  - 


^6 

3  <p 


(2-4) 

(2-5) 


(2-6) 


H  r 
r 


3E 

1  i  1  3  ,-r  \ 

r  — -  Ta  (EAsin(9)  -  rr- ) 


(2-7) 


jwU  rsin9  39  3$ 

respectively.  The  r  component  can  be  neglected  in  many  practical  antenna 
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coupling  situations,  but  we  shall  retain  it  in  the  theory  and  equations  for 
the  sake  of  completeness. 

The  spherical  polarization  components  of  Antenna  A  and  Antenna  B  may 
be  written  as 


0 

cos(e)cos($ ) 

cos(0  )sin(<jj ) 

-sin(0)“| 

X 

* 

- 

-sin($) 

cos ($) 

y 

r 

sin( 9)cos ($) 

sin(0)sin($) 

cos  (<j>)  1 

z 

,  and  (2-8) 

0* 

cos(  0'  )cos(  ) 

cos  (  9*  )sin(  ) 

-sin(  0’) 

x 

♦  ’ 

as 

-sin($ ' ) 

cos ($ 1 ) 

0 

y’ 

sin(0 '  )cos  (d> ' ) 

sin(9  *  )sin(<j> '  ) 

COS  (cj>  '  ) 

i' 

(2-9) 


respectively.  The  '  symbol  is  placed  over  the  near-field  0  and  <j> 
variables  to  distinguish  them  from  far-field  angular  variables  that  will  be 
presented  later.  The  (x,y,z)  and  (x'jy'jZ*)  are  parallel  non-rotating 
cartesian  unit  vectors  fixed  at  the  center  of  rotation  of  Antenna  A  and 
Antenna  B,  respectively.  The  angles  5  and  £  for  Antenna  A  may  be 
expressed  as  functions  of  the  angles  §'  and  for  a  point  on  the  surface 
enclosing  Antenna  B.  In  particular,  the  expressions  for  sin(0)  and  sin(<j>) 
are 


sin( 0) 


sin(  <j>) 


1  - 


Rcos(0  )  +  r' cos (O') 
_ o _ 

|R+r’i 


1/2 


,  and 


(2-10) 


Rsin(  0o)sin(  <p0)  +  r'  sin( 6' )sin( $' ) 
i  R  +  r'  I  sin(9  ) 


(2-11) 


We  note  here  that  the  conditions  that  make  sin(0)  equal  to  zero  also  cause 


£ 

15 

> 

s 

f\ 

3 

v 

£ 

r 

<.■ 

V’ 


V 

V 

'j 

V 


■s 

V' 


S 

.*» 
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l' 

7S 


n 


the  numerator  of  Equation  (2-11)  to  equal  zero.  More  precisely,  the 
maaerator  and  denominator  of  Equation  (2-11)  both  approach  zero  as  sin(6) 
approaches  zero  and  the  proper  limiting  value  of  Equation  (2-11)  is  zero. 

The  projection  of  the  polarization  components  of  Antenna  A  onto  the 
surface  enclosing  Antenna  B  can  be  accomplished  via  the  transformation 
matrix  [T] ,  to  wit: 


T 

00' 

T 

A$0 

T  , 

T  , 

0$’ 

T  „  . 

T 

0r 

V 

Tre>~ 

e; 

T  , 

Ea 

rr’ 

r 

(2-12) 


The  elements  of  [T ]  are  computed  by  forming  the  scalar  product  of  the 

A  |  A  A  A 

subscripted  unit  vectors.  For  example,  Tgg»  ■  9  *  9' ,  T^r»  *»  .  r', 

and  so  forth.  The  elements  are  tabulated  in  Table  2-1  for  convenience. 

The  vector  operations  contained  in  the  integrand  of  the  coupling 
integral  can  now  be  performed  with  th»  aid  of  Equations  (2-6)  through  (2- 
12).  Let  the  surface  enclosing  the  receiving  antenna  be  a  sphere  of  radius 
r'.  Thus,  the  unit  normal  n'  to  the  surface  is  then  identical  with  r’ . 
Consequently,  only  the  r*  component  of  the  vector  cross  product  terms 
within  the  square  brackets  of  the  integrand  survive  when  the  scalar  product 
is  formed  with  r' .  The  coupling  integral  is  written  as  the  sun  of  four 
integrals, 


,ba 


-«E 


HbJ> 


«Ea 


H>  - 


((H 


8’ 


♦  ((H* 


E> 


(2-13) 


where  the  integrals  are 


((Ee”HJ'»-y/lEe  a*!’>  ree’  *  <***'>  V 

S’ 

♦  E*(R+?’)  TrQl]Hb,(F’)  ds', 


(2-14) 
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MATHEMATICAL  EXPRESSIONS  FOR  THE  ELEMENTS  OF  THE  POLARIZATION 


/ 


/ 


«E*"H8'»  -  fM  (5"'>  *  E‘,  (5*''>  V 

S’ 

♦  E*(R+r ' )  Tr^,)Hg,(r')  da', 


(( 


J  ■'c* 


[H“  (R+T* )  T00,  +  Ha  (R+f ' )  T^ 


+  H*(R+f ' )  da',  and 


■//. 


[Ha0  (R+?' )  T04 
S’ 

*  h£(R*?')  T^J 


,  ♦  H*  (R+f) 

)  $ 

Eg,(f )  da’  . 


(2-15) 


(2-16) 


(2-17) 


These  integrals  express  the  coupling  in  terms  of  the  spherical  polarization 
components  of  the  near- field  E  and  H  fields  of  the  two  antennas  and  the 
components  of  a  spherical  polarization  transformation  matrix.  Similar 
expressions  can  be  derived  for  cartesian  polarization  components,  but  they 
are  not  needed  in  this  analysis. 

C.  SAF  Coupling 

The  integrals  (2-14)  through  (2-16)  require  the  use  of  the  near-field 
E  and  H  fields  of  the  two  antennas  to  compute  coupling.  The  near-field  E- 
and  H-fields  for  most  antennas  are  not  ordinarily  known.  They  can,  of 
course,  be  measured  or  they  can  be  computed  from  the  known  far-field 
quantities.  Thus,  considerable  time  and  money  resources  would  have  to  be 
devoted  to  obtain  the  near-field  data.  Then  one  still  has  to  perform  the 
integrations  (2-14)  through  (2-17).  Thus,  it  will  be  advantageous  to 
express  the  integrals  in  terms  of  the  far-field  E  and  i?  fields  of  the 
antenna . 

Each  of  the  four  coupling  integrals  (2-14)  through  (2-17)  is  the  sum 
of  three  coupling  integrals,  all  of  which  have  the  same  general  form.  In 
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a  W 

particular,  the  coupling  integral  ((EQf,Hfi,))  may  be  written  aa 


((Ee»  'V»  ’  «Ee  Tee'  ’<’»  +  ((Ee“  V ’V  +  ((E!  Tre'  *<-  }) 

Attention  ia  focuaed  on  ((  E*  Tfl ,H^, ))  which  ia  typical  of  the  integral8 

0  00  $ 

under  conaideration.  This  integral  ia 


«  ‘l  Te?-V» 


•  yy^f)  nj.<f|!Te 


„  ds', 


(2-18) 


Since  the  surface  has  been  chosen  to  be  a  sphere  of  radius  r'  centered  on 
Antenna  B,  this  can  be  written  more  explicitly  as 


«Ee  t99-V)  *  (r'y2  f  f  "J<?' 


)  TQOt  sin(0')d 4>'d0' 
00 

(2-19) 


where  it  should  be  noted  that  Tno,  is  also  a  function  of  r'  and  the 

00 

integration  variables  ( 0' ,$'). 

The  integrals  (2-18)  or  (2-19)  belong  to  a  class  of  generalized 
convolution  integrals  defined  on  an  arbitrary  closed  surface.  The 
generalized  convolution  integral  is 


:±(r)  -  yy 


A(R*r ' )  •  B(F ' )ds 1 


(2-20) 


where  the  notion  of  convolution  has  been  extended  to  include  integrands 
involving  A(R+?' )  .  5(r')  as  well  as  the  usual  A(R-F')  •  B(r').  Integrands 
of  the  type  A(F'-R)  •  B(F')  can  also  be  included,  but  they  are  not  needed 
here.  The  generalized  convolution  integral  C  (R)  has  the  useful  property 


ffiat 


?')  •  B(F')ds'-  -2*2 


?/// 


exp  L+Ik  *R1 

2  u2 

k  -  k 


i(K)*b(K)dric 


(2-21) 
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where 


ms  .  *.//  5T  -  1  e3<'r"  d5">  (2'22) 


// 


B(?' )e+J  K ’  da'. 


(2-23) 


- —  is  the  normal  derivative  on  the  surface  SM  where  SM  is  a  surface 
3n  _  _  — 

enclosing  the  sources  of  A,  S'  is  a  surface  enclosing  the  sources  of  B  ,k 

is  a  wavevector,  and  T  ia  a  three-dimensional  volume  in  wavevector  space. 

The  relation  (2-21)  holds  under  conditions  that  are  much  less 

restrictive  than  those  encountered  in  real-world  antenna  coupling  and 

scattering  problems.  Sufficiency  conditions  for  (2-21)  are  (1)  the  vector 

fields  A  and  B  each  satisfy  the  homogeneous  vector  Helmholtz  equation 

outside  their  respective  source  regions,  (2)  the  sources  of  A  are  contained 

inside  a  finite  surface  s  <  S"  and  the  sources  of  B  are  contained  inside  a 

a 

finite  surface  sb<  S',  and  (3)  S'  and  S"  do  not  overlap. 

The  derivation  of  (2-21)  is  readily  accomplished  under  the  conditions 
just  stated.  The  field  A  (R+f)  can  be  expressed  as  the  Huygens-Fresnel 
integral  [17] 


A(R+f ' )  (  A(f")  ~  - 


3Akp  ]  G(R+T '  |  f")  d?", 


(2-24) 


where  denotes  the  normal  derivative  at  the  surface  S"  and  where  G  is 

the  free  space  Green's  function  given  as 


G(R+f '  |  ?»)  ,  exp (- jk  |(R+? ' )  -  j  ] 
j (R+r ' )-r" j 


(2-25) 


/ 


Next,  Che  coordinate  space  represenCaCion  of  G  is  replaced  by  Che 
wavevector  space,  or  spectral,  representation  Gk  given  as 


Gk(R  ♦  f '|  ?"  ) 


*/// 


*  [-jtc.(R+r'-r"  )] 

e*P _ _  dr 

2,2  tc 

k  -k 


(2-26  ) 


where  f  is  a  three-dimensional  volume  in  wavevector  space [18  ].  The 
"shape”  of  f  and  the  integration  limits  will  be  specified  shortly.  If  Gk 
is  substituted  into  (2-24),  and  (2-24)  is  substituted  into  (2-20),  we 
arrive  at  the  expression 


C*(5) 


f 

hr  f/j f  h  -  ^ 

r 

t-  e" 


]  exp[j<*r"]  ds" 


If 


B(r')  exp[+jK*r']  ds'  >  dT^  , 


(2-27) 


which  is  equal  to  the  right-hand  side  of  (2-21).  The  order  of  the 
integration  was  interchanged  in  arriving  at  (2-27).  This  is  permissible 
because  Che  integrands  involving  A  and  B  are  "square  integrable"  over  their 
respective  domains.  The  integrands  involving  A  and  B  are  square  integrable 
because  they  are  produced  by  sources  of  finite  extent.  For  the  integral 
over  a  cont°ur  which  circumvents  the  simple  poles  of  the  wavevector 
integral  located  at  k  m  *k  is  chosen  in  order  to  properly  represent  the 
Green's  function  (2-25).  The  wavevector  integrand  is  finite  and  integrable 
along  the  specified  path. 

It  will  be  advantageous  to  use  the  spherical  integral  representation 
of  the  Green's  function  given  as  [18,19] 


•  %  % 
v.v* 

>\%\N 


V*./ 

V/./. 

W* 

»v_y 


G(R+rM  r") 


Iff 


exp[-j<’  (R+?'-r")]  2 

_ _  K  SI 


2  i2 

k  -k 


neditd  0d< 


(2-28) 
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ft*  -i*  *•**»*  ** *  ***•-* 


i -*•*  nil •y-**»‘-n •»•*»  '* 


*_*  '»  •  *  nn  * k  .■'i  r  >  **  k* 


va  •a.^-4'^-ji •‘*J»  V  V 

& 

■f.V 


where  the  integration  over  k  follows  the  contours  shown  in  Figure  2-2. 

Thus,  the  path  of  integration  is  along  contour  Cj  for  cosfy)*  0  and  along 
countour  C2  for  cos(y)  >  0,  where  y  denotes  the  angle  between  <  and  R. 

Since  functions  a(ic)  and  b(x'  are  analytic,  the  only  poles  are  the  ones  at 
k  *  *k.  Hence,  application  of  the  residue  calculus  to  Equation  (2-27}  yields 


IT  IT 

:*(R)  -  k  //  exp[-jic(9,$)  *  R  ]  a  (6,$)  *  b(0  ,$)sin0d<|>d0 


(2-29) 


where 


k(0,if>)  *  kosin(0)cos($)  x  ♦  k0  sin(9)sin(<fr)  y  ♦  k0  cos0z  (2-30) 

This  result  is,  of  course,  applicable  to  each  scalar  component  a  •  b  of  the 
integrand. 

The  partial  coupling  integral  of  Equation  (2-19)  can  be  expressed  with 
the  aid  of  Equation  (2-29)  as 


«  E0a  T00-  >HJ<» 


7T  TT 

If 


exp[-jk(0,$)  .  R]  F®(0  ,<p  > 


IT,  (9-8  ,  <p— a  )sin(0)  d<t>d9 


(2-31 } 


where  8  and  a  are  the  rotation  angles  depicted  in  Figure  2-3, 


F? (9  ,<p  )  -  r  eJkr  (0,*),  and 


(2-32) 


,b  , 


where  U^,  is  a  functional  of  the  H  , -field  of  Antenna  B  given  by 


(r  |2/  /  t"*’  <r’’5’ 


-8,<J>'-a) 


'  (1  '  ,e' ,$' )  exp[- jk(e  ,<p )  .  f]  sin(0' )|dVd9' .  (2-33) 


U^r (9-6 ,0-®) 


The  generalized  convolution  integrals  discussed  previously  can  be  utilized 

along  with  an  approximate  relation  applicable  for  E  «  (E. ,  Ex), 

r  0  9 


(2-34) 


where  Zc  is  the  impedance  of  free  space, to  express  (2-33)  as  U^,  -  U„,  , 
b  <f  a 

where  U  ,  is 

0 


Ug,  (0-6  ) 


k  b 

■yf  F0t  (5-6,n-  a  )  T06(  (0 , j  C,n)sin£dnd£ 


o  -ir 


(2-35) 


The  angles  (£  ,n )  correspond  to  the  usual  (0  ,$  )  spherical  angles.  The 

functions  F**,  and  T  are 
0  00 


Fg,  (€-S,n-a  )  -  rejkr  Eg,(C-B,n-a)  and 


(2-36) 


IT  IT 

Te0^e  >$U  »n  )  *  (r')2  y  y*  exp[-j[k(0,<(>)+  k(c,n).l*  ]  TQQ,  dsi  ’ 

o  -ir  (2-37) 


respectively,  where  d«'  =■  sin(S' )d§ 'd$' . 

Equation  (2-31)  can  now  be  written  more  explicitly  as 


«  Ee  Tee  .-»*-» 


2  7T 


If 


exp[-jk(0,4>).  R]  FQ ( 0 , d>)  sin(0)  x 
0 


//  F0»  (5-B,n-a  )  T09,  (9,4>|€,n)s inC  dnd£  d<j>d9 


(2-38) 


-  ’v v  v.* v  v v  v  v  v  v v.* 


This  expression  for  the  partial  coupling  integral  has  the  receiving  antenna 
SAF  Fb,  convolved  with  the  integral  transform  TQq,  of  the  polarization 
projection  function  T0g,.  Thus,  this  convolution  must  be  performed  for 
every  new  rotation  angle.  However,  for  our  basic  computation  scenario  we 
wish  to  rotate  the  receiving  antenna  while  keeping  the  transmitting  antenna 
stationary  during  this  set  of  computations.  Accordingly,  it  would  be 
computationally  more  efficient  to  have  an  expression  in  which  the  non¬ 
rotating  SAF  Fa  of  the  transmitting  antenna  rather  than  Fb  of  the  receiving 
9  w 

antenna  was  convolved  with  Tggl.  This  is  readily  achieved  simply  by 

interchanging  the  order  of  integrations  in  (2-38).  The  dunmy  variables 
( 9,  <p)  and  (  £,  n)  are  also  exchanged  after  the  re-ordering  of  the  integrals  so 
that  the  final  integration  is  expressed  in  terms  of  the  familiar  (0,0) 
spherical  angles.  Thus,  the  coupling  integral  becomes 

r  j  j  Fb0,  (9-0,0  -a)  sin(0)  x 

O  -IT 

ir  it 

F6(C’ri)T00'(0^l5’n>exp["^(?*Tl)* 

O  “IT 

sin(C)  dnd£d0d9  .  (2-39) 

*  e  K 

The  equivalence  of  (2-381  and  (2-39)  is  expected  since  ((  EQTQQ,  ,H,))  is 

O  DO  9 

just  a  multiple  convolution.  Equation  (2-38)  is 


«  Ee*V  •  »]  *  rj  *  V  (2-40) 

where  the  *  here  denotes  convolution.  Similarly,  Equation  (2-39)  is 

«  EeTee'  ’$■  »  ‘  f6  *  C«pf-Jk.Slr*)  *  t96, 

Sufficiency  conditions  for 
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(2-42) 


P*  *  (exp[-jk*E]Pg)  *  T0QI  -  (exp [- jk *R]F®)  *  ?*  *  T00, 


are  Che  saae  as  given  previously  in  Che  discussion  of  Che  generalized 
convoluCion.  Less  rescriccive  condi Cions  for  Che  commutability  of 
convolution  integrals  are  given  by  L.  Schwartz  in  References  20  and  21. 

The  sane  procedure  is  applied  Co  all  of  Che  partial  coupling  integrals 
except  Chose  having  and  Hr  in  their  incegrand.  These  integrals  are  of 

A  A 

Che  order  r  whereas  Che  integrals  involving  only  produces  of  Che  (  9,  <f>) 
and  (  ©,<#»>  components  are  of  the  order  r  .  The  integrals  involving  E 
and  are  expected  to  yield  a  non- negligible  contribution  to  the  coupling 
only  if  one  (or  both),  of  the  antennas  is  (are)  electrically  very  small  and 
the  separation  distance  is  about  two  wavelengths  or  smaller.  The  emphasis 
in  this  research  program  is  placed  on  development  of  efficient  algorithms 
for  directive  antennas.  It  appears  that  the  contributions  of  the  radially 
polarized  components  will  be  negligible  for  realistic  situations  involving 
directive  antennas  and  for  many  situations  where  one  or  more  of  the 
antennas  are  electrically  small. 

The  total  coupling  integral  ((b,a))  *  V^a  is  then  expressed  as  the 

sum  of  the  double  convolution  integrals  for  all  of  the  partial  couplings 

involving  the  (9 )  components  of  the  transmitting  antenna.  The  expression 
ba 

for  V  is  then 

,b* '  ~k  [■ pb  *  "’■p1-*51"'!’  *  <v +  v> 

♦  Pj  *  (exp(-jk*R] F®)  *  (T99,  +  T99,) 

*  Fe  *  <-P(-iE-5lPj)  *  (T99.  -  t99.) 

-  pj  .  (exp[-jk*R]Fg)  •  (!„,’+  T99,)J 

(2-43) 
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The  antenna  coupling  ia  more  conveniently  characterized  in  terms  of  the 
mutual  gain  M  for  the  antenna  pair  relative  to  a  pair  of  isotropic 

Ka 

radiators  in  lieu  of  the  voltage  V  .  The  mutual  gain  ia  defined  as 


M(R,8,a)  - 


Pr(R,  6,o) 


(2-44) 


where  Pr  is  the  received  power  and  is  the  transmitted  power.  The  mutual 
gain  is  expressed  in  terms  of  the  normalized,  dimensionless  SAP's  Fq  and  Fq 
that  are  normalized  to  unity  as 


H(R,e,o)  -  8b8a  R4  |^0  *  (expI-jk-RlF^)  *(T00,  +  T^,) 

*  C  *  *  (Tee’  +  V* 

*  *  <V  +  V’ 

-  ^  *  (.xPl-jk  •!»;„>  ♦  (tm,  .y  2 

(2-45) 


where  Sb  denotes  the  peak  "far-field"  gain  relative  to  isotropic  gain  of 
Antenna  B  when  used  as  a  receiver  for  plane  wave  illumination  and  ga 
denotes  the  peak  far-field  gain  of  Antenna  A  when  used  as  a  transmitter. 

The  mutual  gain  given  by  this  equation  is  the  quantity  of  interest  in  the 
numerical  studies. 

D.  Obstacle  Effects 

We  consider  next  two  antennas  operating  in  the  presence  of  an 
obstacle.  The  obstacle  is  allowed  to  have  arbitrary  but  known  shape  and 
composition,  and  it  can  be  located  anywhere  outside  the  surfaces  enclosing 
the  radiating  portions  of  the  two  antennas.  The  derivation  of  the  coupling 
integrals  for  two  antennas  operating  in  the  presence  of  a  single  obstacle  is 
outlined  in  the  following  paragraphs.  The  extension  to  more  than  one 
obstacle  is  then  briefly  considered. 
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The  situation  under  study  is  depicted  in  Figure  2-4;  R  is  again  the 
vector  extending  from  the  transmitting  Antenna  A  to  the  receiving  Antenna 
B.  Rg  and  r  are  vectors  extending  from  Antenna  A  and  the  obstacle 
reference  center,  respectively,  to  a  field  point.  ra  and  rb  extend  from 
Antenna  A  and  Antenna  B,  respectively,  to  the  obstacle  reference  center,  r^ 
traces  out  the  surface  of  the  obstacle. 

It  is  well  known  that  the  E  and  H  fields  at  a  point  in  the  near  field 
can  be  written  as  the  sum  of  the  incident  antenna  fields  and  the  obstacle 
scattered  fields, 

Ec  -  Ea  +  Ea,  and  (2-46) 

Hc  -  H«  +  Ha,  (2-47) 


where  the  superscripts  t,  a,  and  s  denote  total,  Antenna  A,  and  scattered, 
consecutively.  The  complex  voltage  Vbt  induced  in  Antenna  B  may  be 
obtained  from  the  coupling  integral  given  in  Equation  (2-3)  by  substituting 
(Ic,  Hc)  for  (Ea,  Ha)  to  yield 


ybt  m  yba  +  ybs 


(2-48) 


The  voltage  Vba  is,  of  course,  just  the  voltage  induced  in  Antenna  B  by 
Antenna  A  in  the  absence  of  the  scattering  obstacles,  i.e.,  the  clear-site 
coupling  voltage.  Thu  voltage  Vbs  is  the  voltage  induced  in  Antenna  B  by 
the  scattered  fields  of  the  obstacle,  and  it  is  given  in  terms  of  the 
coupling  integral  as 


ybs 

1 


(-fb+f')  x  Hb(f ' )  +  Hs  (-rb+?')x  Eb  (r') 


n'  ds'  ] 
(2-49) 


where,  as  before,  S'  is  an  arbitrary  surface  enclosing  Antenna  B  and  h'  is 
a  unit  vector  normal  to  the  surface  and  pointing  out  of  the  enclosed 
volune.  Proceeding  exactly  as  in  the  derivation  of  the  clear-site  SAF 
coupling  integrals,  the  voltage  Vbs  is  expressed  in  terms  of  the 


<-■  v  -/ 


V  vv’ s*  r-T  C  V  *«'  V  O  *.*  \  *.* 
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+  FJ  *  <exp[-jk-R]F^  )  *  (T9QI  +  Tw§) 


+  Fe  *  (exp[-jk-R]F^  )  *  (Tq<>,  +  T^,) 


•  Ha  (ra  +  rpldsi 


(2-51) 


c&  \ 

-J  ■ 

.<2  ; 


far-field  quantities  Fb(®,$)  and  Fa(®,$)  as 


ubs*  — 5- 

1  Z 


F0  *  <e»pt-jk-RJF*e  )  •  (Tm,  +  Tm,) 


-  *  (.xpl-jk-WF’  )  *  <T„,  +  T 


,j 

J(2-50) 


The  SAF  Fj(9,<)>)  for  the  scattering  obstacle  may  be  obtained  with  the 
aid  of  a  slight  generalization  of  the  reaction  integral  for  scattering 
[15].  Let  Jje(?1)  and  )  denote  scattering  current  dyadics  for  the 

electric  and  magnetic  currents  induced  on  the  surface  of  the  obstacle  when 
illuminated  by  an  infinitesimal  dipole  with  unit  current  located  at  ?  * 
(r,0,<j>).  In  the  limit  of  very  large  r,  we  obtain,  via  the  reaction 
integral  for  scattering 


T(  e,  <» 


^  If  [v  •  «v  •  5‘ 


+  r  )-J  (r  ) 
1  1®  1 


from  whence  we  can  obtain  F^C  9,  <J>)  via  the  relation  Fa(  0,  <p)  •  reJ^r  Esj  (  9>  ♦) . 

An  expression  for  F9  (  9,  d>)  in  terms  of  far-field  pattern  variables  is 
desired,  and  can  be  obtained  by  employing  the  Huygens-Fresnel  integral 
given  by  Equation  (2-24).  In  particular,  the  electric  field  Ea  of  Antenna 


A  becomes 


E*(ra  *  v  -  £  If  Fa( 0', 


$')  exp[  -jk(0',<j),)*(?a+?i)  j  dC'  , 

(2-52) 


and  similarly  for  H®.  This  leads  to  the  following  expression  for  F®  (9,#), 

?!<M)  **lrff  fa<0v>exp[-jfc<6' >*’>•**] • 

^^*31(r1)exp  [-Jk(0',*,)*r  J  ds^n* 


(2-53) 


where  J,  *  J,  +  J,  and  where  dH  '  ■  sin(9')  d0'd$'.  This  can  be  cast 
l  le  lm 

into  a  more  familiar  form  by  defining  the  obstacle  scattering  dyadic 
^<0,*!  9',*'), 

3^  (9  >$ |  9  * ,♦  * )  *  j«li  ff  (Fj  1 0  ,$)exp  |  "i^(e  '  »$')•  dsx  . 
1  1  (2-54) 


The  resulting  expression  for  F*1(0,^)  is  then 


F®(0,$)  “  — 
X  T  2irwp 


f f  s  Ce ,4|e*«4*>  •  Fa(9\*')  exp[  -jR(0',4>  •  ?a]dn  • 

—  J  ^  /A 


(2-55) 


This  equation,  as  well  as  Equations  (2-51)  through  (2-54)  were  derived  with 
respect  to  polarization  components  that  are  convenient  for  the  description 
of  the  surface  currents  on  the  obstacle.  However,  Equation  (2-55)  is  valid 
for  any  choice  of  polarization  convention.  Accordingly,  spherical  (9,<M 
components  are  now  chosen  in  order  for  F®(e,$)  to  blend  with  the  "dear- 
site"  SAF  equations. 

The  foregoing  analysis  is  readily  extended  to  multiple  obstacles  if 
the  multiple  obstacle  scattering  dyadic  S  can  be  specified  from  theory  or 

M 

experiment.  S  can  be  approximated  as  the  sun  of  the  isolated  single 
obstacle  scattering  dyadics  if  the  effects  of  multiple  scattering  can  be 
ignored.  Perturbation  techniques  for  small  metal  cylinders  or  spheres 
exist  for  improving  the  results  obtained  from  the  superposition  of  isolated 
obstacle  scattered  fields  [22,23]. 
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SECTION  III 

APPROXIMATE  SAP  ANALYSIS 


A«  Introduction 

The  rigorous  theory  and  equations  presented  in  Section  II  have  the 
very  desirable  attribute  of  being  applicable  to  antennas  and/or  scatterers 
of  low,  moderate,  or  high  directivity.  Additionally,  they  are  valid  for 
any  orientations/ arrangements  of  the  antenna  and  scattering  objects. 
Accordingly,  the  rigorous  SAP  equations  appear  to  be  well  suited  for 
analyzing  antenna  coupling  under  both  clear-site  and  obstructed  conditions 
for  situations  where  there  are  questions  concerning  the  PWS  computations. 
However,  one  can  anticipate  that  the  computer  core  storage  and  execution 
times  for  electrically  large  antennas/sc.,  terers  may  become  prohibitively 
large  for  a  computer  program  based  on  the  rigorous  SAP  analysis. 

The  rigorous  SAP  antenna  coupling  analysis  requires  the  evaluation  of 
complex  integral  transforms  of  (T„„  ,  T, .  ,  etc.)  of  the  polarization 

00  (p<p 

transformation  operators  (T„„ ,  T, ,,  etc.)  and  the  evaluation  of  a  pair  of 

00 

nested  convolution  integrals.  The  integral  transform  is  a  four-dimensional 
complex  function  of  angular  variables  while  each  of  the  convolution 
integrals  are  two-dimensional  complex  functions  of  angular  variables . 
Assuming  that  the  functions  are  sampled  at  the  Nyquist  rate  over,  say,  a 
hemisphere,  it  can  readily  be  deduced  that  computer  storage  and  execution 
times  can  become  excessive  for  closely-spaced  electrically-large  antennas. 
Indeed,  storage  requirements  and  execution  times  are  significant  problem 
areas  even  for  the  inherently  fast  PWS-based  algorithms  that  employ  FFT 
processing  [10,11].  Accordingly,  it  is  recognized  forthwith  that  some 
simplifications  must  be  made  in  order  to  obtain  a  practical  and  versatile 
SAP  algorithm  that  can  be  applied  to  a  large  class  of  antenna  coupling 
situations . 

Analytical  and  numerical  studies  were  initiated  during  this  reporting 
period  to  improve  the  efficiency  of  the  computer  programs  employing  SAF 
equations  without  seriously  degrading  accuracy.  The  stationary  phase 
evaluation  of  radiation  integrals  has  been  applied  along  with  sampling 
theory  to  furnish  guidelines  for  choosing  the  required  angular  integration 
regions  for  the  convolution  integrals  and  the  integral  transforms.  The 
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stationary  phase  analysis  was  also  employed  in  the  development  of  series 
representation  of  the  integral  transforms  TQ0,  T ,  etc.  Numerical  studies 
were  begun  to  study  the  impact  of  integration  region  and  sample  spacing  on 
the  accuracy  of  the  antenna  coupling  computations.  The  analysis  and 
preliminary  numerical  data  are  presented  and  discussed  in  the  following 
paragraphs . 


Stationary  Phase  Evaluation 


ration  Region 


The  integral  transforms  such  as  TQ0  given  in  Equation  (2-37)  and 
the  convolution  integrals  such  as  ((E|  TQQ,  H^,))  given  in  Equation  (2-39) 
are  of  the  general  form, 


I<Yi,Y2)  *  ff  A(Y1»YP  «*p[“jO«)  *  <Yi»Y2|y{,  Y2)j  dY^dY* 

(3-1) 


P.  Debeye  [  24  ]  has  devised  a  method  for  evaluating  this  type  of  integral, 
known  as  the  method  of  stationary  phase,  or  the  method  of  steepest  descent. 
The  method  is  based  on  the  observation  that,  for  kr  >>1,  the  major 
contribution  from  the  integrand  comes  from  the  region  in  (yJjYp  space 
where  the  phase  function  (kr)  f ( Y j , Y2 1 y[,Y2)  i»  slowly  varying  compared  to 
the  largest  significant  variation  of  the  function  A(Yj,Y?).  This  region  of 
relatively  slow  phase  variation  is  centered  at  the  point  where  the  first 
derivative  of  the  phase  function  vanishes,  known  as  the  stationary  phase 
point.  The  region  of  comparatively  slow  phase  variation  is  frequently 
referred  to  as  the  stationary  phase  region. 

The  integrand  rapidly  becomes  increasingly  oscillatory  when  the 
integration  extends  outside  the  stationary  phase  region.  The  contribution 
of  this  exterior  region  is  typically  very  small  compared  to  the 
contribution  from  the  stationary  phase  region,  provided  the  maximum  of  a 
sharply  peaked  A(Y[,Y2)  i3  located  well  within  the  stationary  phase  region 
or  provided  that  the  parameter  kr  is  sufficiently  large. 

In  order  to  find  the  stationary  phase  point,  the  phase  function  kr 
.  »  » 

f (Yj, Y2 1 Yj,Y2)  is  expanded  to  second  order  in  a  Taylor  series  in  two 

»  » 

variables  about  the  points  Y^o  and  Y20, 


kr  •  f<Yi,  Y2 1  Yi  o  +h,Y2o  +  g)  -  kr  •f(Yl0*Y20>  + 

h~*s  -Ar]  (kr)  f  <  Yi,  Ya  i  Yi  ,  Y2)  |  ’  ' 

3Yi  3Y2J  10  20 

2 

h  -K  *  g  -V]  f  (kr)  f(Y1,Y2|Y[  ,Y2)  I  ,  ,  , 

3Yj  3Y2  J  Y10  *  Y20 

(3-2) 

and  Che  stationary  phase  point  is  found  from  the  requirement  that 

fh  -A-  1  g  k0r  f(  Yi,  Y2|  y[  ,  Y2)  |  *  *  (3-3) 

L  3Yi  3Y2J  Y10*  Y20 

The  phase  functions  to  be  evaluated  in  the  SAF  integrals  can  be  written  in 
the  form 

(kr)  f(Yi,Y2|Yl  ,  Y2)  a  (kr)  {sin(  Yj  )cos  (  Y2)sin(  yJ  )Cos  (  Y2) 

+  sin(  Yi)sin(Y2)sin(Y1)sin(Y2) 

+  cos(Yj)cos(Yi)  }  ,  (3-4) 

Expanding  the  phase  function  in  its  Taylor  series  expansion  and  applying 
the  condition  (3-3)  leads  to  the  sensible  result  that  the  stationary  phase 
point  is  located  at  Yj  *  Yj  and  Y 2  »  Y  2  . 

The  extent  of  the  stationary  phase  region  can  be  estimated  as  follows  [11]. 
First,  a  coordinate  rotation  is  effected  to  align  the  z-axis  of  a  spherical 
coordinate  system  along  the  stationary  phase  direction  (Yj,Y2>.  The  phase 
function  in  this  new  coordinate  system  is  then 

(kr)  f ( •©■ )  =  kr  cos('90  (3-5) 

where  &  is  the  polar  angle  that  opens  from  the  new  z-axis.  Next,  the 
criterion  is  adopted  that  the  edge  of  the  stationary  phase  region  occurs  at 
the  angular  location  where  the  rate  of  change  of  the  phase  function  is 
equal  to  twice  the  most  rapid  variation  of  the  function  A.  Thus,  we  need 
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Co  equate  Che  firsC  derivacive  of  Che  phase  function  with  respecc  Co  Che 
angular  incegracion  variables  Co  Cwice  Che  highesC-frequency  Fourier 
componenc  of  A.  For  (non-supergain)  anCenna  paCCern  functions  and/or 
obaCacle  scaCCering  functions,  Cwo  times  Che  highesC-frequency  Fourier 
componenc  of  A^)  is  bounded  by  tdA  ,  where  D  is  either  (1)  Che  diameter 
of  che  sphere  circumscribing  Che  radiacing  portion  of  Che  anCenna  or 
scacterer  or  (2'  ♦■wo  wavelengths,  whichever  is  larger.  For  coupling 
integrals,  for  which  AC^1)  is  essentially  a  product  of  two  non-supe.'gain 
anCenna-like  pattern  funcCions,  Cwice  the  highesC-frequency  component  can 
be  taken  as  w  (Da  +  Db)/A  where  Da  and  Db  are  (1)  either  Che  diameter  of 
spheres  circ unscribing  the  radiacing  portions  of  Antenna  A  and  Antenna  B, 
respectively  or  (2)  two  wavelengths  each,  whichever  is  larger.  Thus,  for 
coupling  integrals  we  have 

(kQr)  sinOO  -  ir  JMRbl,  (3-6) 


where  -Q^  denotes  the  "edge"  of  the  stationary  phase  region.  This  leads  to 
the  result  [11] , 


sin(  &  ) 
P 


(Da+Db) 

2r 


(3-7) 


The  sample  spacing  can  be  chosen  by  invoking  the  Nyquist  criterion. 
Accordingly,  the  angular  increments  should  not  be  larger  than 


a  y  = 


l 


IT  fD  +D  ) 
a  b 


(3-8) 


in  order  to  satisfy  the  Nyquist  sampling  rate. 

The  integ-  it  on  region  and  sanple  spacing  given  by  Equations  (3-7)  and 
(3-8),  respectively,  provide  guidelines  for  numerical  computations.  Larger 
integration  regions  and/or  finer  sample  spacings  may  be  required  in  some 
applications  in  order  to  achieve  a  specified  level  of  accuracy.  On  the 
other  hand,  some  relaxation  of  these  criteria  may  be  permissible  for  some 
EMC  applications  The  numerical  studies  that  are  in  progress  on  this 


contract  should  provide  additional  insight  into  the  impact  of  the 
integration  region  and  sampling  rate  for  EMC  applications. 


2.  Approximation  for  the  Integral  Transforms 

An  approximation  is  sought  for  the  angular  space  representation, 
i.e.,  the  integral  transforms  T00,,  T^„  etc.,  of  the  polarization 
transformation  operator.  It  will  suffice  to  consider  only  the  integral  for 
T00 i since  the  same  procedure  is  applicable  to  all  of  the  T's.  The  integral 
for  Tg0  is  given  by  Equation  (2-37)  as 

IT  IT 


T00(e,*U,n)  -  (r 


>2// 


exp[-j  E(0,$)  +[K(S,n)]*  r']  T 


00 


0  -IT 


where  T  *  is  given  in  Table  2-1  as 
00 


sinCe’Jd^'dO 

(3-9) 


T00'  *  cos (0)cos (5)coa (0 ' )cos (3 ' )  ♦  cos(0)sind)cos(§')sin($’) 


♦  sin(0)sin(0’) 


(3-10) 


The  integrations  in  Equation  (3-9)  are  performed  over  the  surface  of  a 

sphere  circumscribing  Antenna  B.  Axonometric  plots  of  Tgg'  and  T^'  are 

shown  in  Figures  3-1  and  3-2,  respectively,  for  a  circumscribing  sphere  of 

radius  5A  ,  where  X  is  the  wavelength.  These  plots  display  the  functional 

form  of  T00*  and  T^>  versus  elevation  angle  0'  and  azimuth  angle  3' . 

In  order  to  develop  a  conventional  asymptotic  series  representation 

for  T  i,  it  is  necessary  to  expand  the  exact  polarization  transformation 
00 

operator  T  •  in  a  two-dimensional  Taylor  series  about  the  stationary  phase 
point  of  the  phase  function  [E(  C,  n)  +  E(0,  $)]•  ?(§',$').  The  exact  phase 
function  is  also  replaced  by  its  two-dimensional  Taylor  series  expansion  to 
second  order  about  the  stationary  phase  point,  which  occurs  along  the 
direction  of  the  resultant  wavevector  K  *  E(  5,  h)  +  E(6,  $) .  This  leads  to 

an  approximation  for  T  •  of  the  form 

00 


^001  ^  5,  n  |  9,  <p)  M  exp  {  - j kr'  [dj^xo  +  ^y^yo  +  ^z^zo^  ^ 
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Axonometric  plot  of  the  polarization  transformation  function 
for  a  sphere  of  radius  equal  to  5  wavelengths. 


where  are  the  coefficients  of  the  power  series  expansion  for  T  • 
sin(§')  and  where  bj  and  b2  are  computed  as 


«>1  "  -~7  t  Mi  ♦  UyUy  ♦  UZUZ]  |g,  i,  ,  and 
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-  f  -  t  , 


b2  "  — ,  I  ux«x  ♦  Vy  +  °2U^  1 6'  5* 


(3-12) 

(3-13) 
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The  variables  Ux,  Uy,  Uz,  are 


Ux  -  sin(£)cos(n)  ♦  sin( 0)cos ( $) , 


(3-15) 


Uy  ■  sin(£)sin(n)  +  sin( 9)sin( <j>) ,  and  (3-16) 

Uz  ■  cos(0  ♦  cos (6)  .  (3-17) 

Similarly,  the  Ux,  Uy  and  Uz  are 

u;  -  sin<fl')cos($'),  (3-18) 

Uy  ■  sin(9 ' )sin($ ' ) ,  and  (3-19) 

o;  -  cos (S')  .  (3-20) 


The  variables  Uxo,  UyQ  and  Uz0  are  obtained  from  these  last  three  equations 
by  evaluating  them  at  the  stationary  phase  point  (6o,$o). 

The  integral  appearing  in  Equation  (3-11)  is  of  a  well  known  Gaussian 
form.  Closed-form  expressions  can  be  obtained  for  the  various  terms  in  the 
series.  However,  there  are  important  drawbacks  to  this  approximation.  In 
particular,  the  second-order  expansion  for  the  phase  term  is  accurate  only 
over  a  very  small  portion  of  the  stationary  phase  region  unless  kr'  is 
extremely  large.  Many  additional  terms  would  be  needed  in  the  power  series 
expansion  of  the  phao=  function  in  order  to  obtain  an  accurate  approximation 
for  the  phase  function  over  the  stationary  phase  region  for  moderate  or 
small  values  of  kr.  The  integral  would  then  have  to  be  computed 
numerically,  which  would  make  the  approximation  impractical  to  use. 


Moreover,  this  type  of  asymptotic  series  for  T  »  may  not  converge  if  too 

0“ 

many  terms  are  used  in  the  expansion  for  T»  sin(9'),  and  hence  one  must 

00 

determine  the  optimum  number  of  terms  to  use  [25]. 

The  shortcomings  of  the  conventional  asymptotic  expansion  in  the 
variables  can  be  circumvented  by  evaluating  T  *  in  terms  of  the 

f  f  ^ i 

(Ox,0y,Uz)  variables.  Equation  (3-9)  is  thus  rewritten  as 


IT  IT 

Tee’’  <r‘)2  f  f  exp  1  t  Mi  ♦  UyUy  +  u^] }  T00;(u',u')dn'  , 
JJ  (3-21) 
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where 


d  ft '  *  sin(i,)d0,d^'  . 


(3-22) 


The  exact  polarization  transform  operator  T  can  be  written  in  terms  of 
„  00 

(Ui.°z)  « 


T00 '  -[  bA(0i,0i)  u;2  +  B(UX,UX) ]  , 


(3-23) 


where  the  functions  A(Ui,U£)  and  B(U^,U^)  are 

A(Ui,Ui)  -  [( 1+a  Ux)2  -  b  Uz2(  1+a  Ux)  ]  **[  l-U^2^,  and  (3-24) 


B(0]c,Ui)  -[  1+a  0^  -  b  U’2]*  [1+a  S' ]  “*5  [1-Ul2  ]**  .  (3-25) 


The  parameters  b  and  a  are 


,  and 


(3-26) 


R2+r’2 


(3-27) 


The  functions  A(U^,Ug)  and  B(Ui,Uz)  are  expanded  in  a  two-dimensional 
Taylor  series  expansion  about  the  stationary  phase  point  (Uio,U^0)  to  yield 


♦*. 

WaV'.'.C'.V. 


•>  s  s 


a  power  aeries  expansion  for  T  aa»  . 
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^mn  (^x~5jxo^m  (Uz-U20)n 


(3-28) 


Thus,  it  is  seen  from  Equations  (3-28)  and  (3-21)  that  it  is  necessary 
to  evaluate  integrals  of  the  type 

IT  IT 

Imn  -  (  [  Uin  Uin  exp  {-jkr'  { Mx+Vy+0*6z]  }  dA  • 

J  J  (3-29) 

O  -IT 

It  may  readily  be  discerned  that  the  integral  IBn  can  be  obtained  from  the 
integral  I00  for  m  and  n  both  equal  to  zero.  In  particular,  the  integral 
is  obtained  from  I00  by  differentiating  with  respect  to  Ux  and  Uz,  to 


[-jkr  ] 


(nrfn) 


Kl  15] 


(3-30) 


Next,  we  note  that  the  integral  IOQ  obtained  from  Equation  (3-29)  by 
setting  m*n-0  can,  after  a  coordinate  rotation  designed  to  align  the  polar 
axis  in  the  stationary  phase  direction  along  the  resultant  wavevector,  be 
integrated  to  yield 


WUx»uy.uz>  -  (4*) 


sin[kor’U] 

k  r'U 
o 


where 


U  »  [  U2  +  U2  ♦  U2  ]  ** 
x  y  z 


(3-31) 


(3-32) 


and  where  (Ux,Uy,Uz)  are  as  previously  defined  by  Equations  (3-15),  (3-16), 
and  (3-17).  I0Q  is  readily  differentiated  any  number  of  times  with  respect 
to  Ux,Uy  or  uz*  The  derivatives  of  all  order  are  continuous  and  bounded. 

The  general  properties  of  the  (UX,UX)  series  expansion  for  T^i  are 
still  under  study.  However,  preliminary  results  of  the  analysis  indicate 
that  the  series  converges  for  all  r*  <  R,  a  condition  that  is  always 
fulfilled  for  antennas  whose  circumscribing  spheres  do  not  touch  one 
another.  The  preliminary  results  also  strongly  indicate  that  this  series 
converges  much  faster  than  the  conventional  asymptotic  series,  and 
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consequently,  fewer  terns  are  needed  to  achieve  a  specified  acccuracy. 
Work  has  been  performed  to  evaluate  the  CgQ  coefficients  for  the  T  > 
Taylor  series  and  to  obtain  the  analytical  expressions  for  the  integrated 
Iga  terms  out  to  fourth  order  in  0*  and  to  second  order  in  U^.  However, 
the  current  version  of  GTSAF  employs  only  the  leading  integrated  term  I00. 
This  simplification  permits  considerable  reduction  in  computer  run  time 
while  maintaining  good  accuracy. 
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SECTION  IV 
NUMERICAL  RESULTS 


A.  Introduction 

In  Chit  section,  numerical  data  are  presented  and  discussed  to 
illustrate  the  general  capabilities  of  the  SAF-based  antenna  coupling 
algorithm  and  two  related  auxiliary  algorithms.  Firstly,  coupling  data 
were  computed  via  the  computer  program  SSAF  for  a  pair  of  %-A  dipoles  and 
for  a  pair  of  10- X  circular  aperture  antennas.  This  simplified  version  of 
the  full  SAF  computer  program  GTSAF  employs  function  statements  to  generate 
the  spherical  angular  functions  for  the  antennas  under  study.  Also,  the 
full  expression  for  the  integral  transform  Tg0  •  and  i  of  the 

polarization  projection  operator  is  approximated  by  the  leading  term  I0Q, 
given  by  Equation  (3-31),  whose  3-dB  beamwidth  has  been  adjusted  to  more 
closely  replicate  the  mainbeam  of  the  exact  (TqQ  »  +  • ) .  Secondly, 

mutual  coupling  data  is  calculated  via  GTSAF  for  a  pair  of  directive 
antennas  for  clear  site  as  well  as  blocked  conditions  in  the  near-field, 
and  compared  with  experimental  data.  Thirdly,  the  auxiliary  programs  GTOBS 
and  EHPARAB,  for  computing  composite  antenna/obstacle  SAF'a  and  for 
generating  paraboidal  reflector  antenna  SAF's,  respectively,  are  briefly 
discussed. 

B.  General  Trends 

SSAF  was  used  primarily  to  study  general  trends  in  the  mutual  gain 
curves  and  to  study  the  sensitivity  of  the  mutual  gain  curves  to  the  size 
of  the  integration  region  used  in  the  computations.  The  key  results  for 
the  geometry  depicted  in  Figure  4-1  can  be  summarized  with  the  aid  of 
coupled  power  curves  for  a  pair  of  dipoles  and  for  a  pair  of  10-A 

uniformly-illuminated  circular  aperture  antennas.  The  coupled  power 
curves  were  computed  versus  the  elevation  rotation  angle  for  a  fixed 
skew  angle  £0  for  the  geometry  shown  in  Figure  4-1.  It  should  be  noted 
that  the  coupled  power  curves  obtained  for  this  geometry  are  equivalent  to 
the  curves  that  would  be  obtained  when  £0  *  0°  with  receiving  Antenna  B 
located  at  the  longitudal  distance  X  »  R  cos(£0)  and  vertical  separation  Z 
■  R  sin(£0)  and  where  the  rotation  angles  S  and  S',  respectively,  for  the 
two  situations  are  related  as  8  *  S *  —C0 . 
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The  integration  integral  was  always  chosen  to  extend  only  over  the 
uainbeam  of  the  function  (Tgg*  ♦T<^«  ).  The  integration  range  of  the 
outermost  convolution  integral  was  varied  in  the  numerical  tests.  Good 
agreement  was  achieved  between  the  near-field  electric  fields  computed  from 
SSAF  with  one  antenna  replaced  by  an  infinitesimal  dipole  and  the  near¬ 
field  electric  fields  computed  from  aperture  integration  [4]  when 
Equation  (2-38)  was  integrated  in  the  aforementioned  manner.  The 
agreement  was  less  satisfactory  when  Equation  (2-38)  was  integrated  in  this 
approximate  fashion.  Accordingly,  the  SSAF  computations  contained  herein 
were  based  on  Equation  (2-38)  with  the  roles  of  Antenna  A  and  Antenna  B 
reversed.  Since  the  antennas  are  both  identical  and  reciprocal,  the 
coupling  data  computed  in  this  manner  is  equivalent  to  the  coupling  data 
obtained  for  the  situation  depicted  in  Figure  3-3  where  Antenna  A  is  fixed 
at  the  skew  angle  -£0  and  Antenna  B  is  rotated  through  angle  8  .  The  data 
and  discussions  in  the  following  paragraphs  are  presented  in  terms  of  this 
latter  situation  just  described. 

The  coupled  power  curves  computed  from  SSAF  for  a  pair  of  hr  A  dipoles 
separated  by  1A  are  compared  with  the  coupled  power  curves  computed  from 
the  Ohio  State  Piecewise  Sinusoidal  (OHSPS)  wire  program  in  Figures  4-2  and 
4-3.  The  curve  shown  in  Figure  4-2  was  computed  for  a  pair  of  non-skewed 
dipoles,  while  the  curve  shown  in  Figure  4-3  was  computed  for  a  skew  angle 
of  Co  *  -90°.  The  integration  range  for  both  curves  was  the  entire  front 
hemisphere  defined  by  (0°^  9  180°,  -180°^  $  £  180°).  Good  agreement  is 

achieved  between  the  SSAF  and  the  OHSPS  data  for  the  non-skewed 
arrangement.  The  relatively  poor  agreement  for  Co  "  "90  degrees  may 
indicate  that  the  radially  polarized  field  cannot  be  ignored  for  this 
arrangement. 

The  plots  of  coupled  power  shown  in  Figures  4-4  through  4-12  for  two 
10- A  aperture  antennas  provide  useful  and  interesting  information 
concerning  general  trends  and  the  impact  of  the  integration  range  on 
computational  accuracy.  The  figures  are  sequenced  according  to  increasing 
separation  distances  of  10  A  ,  25  A  ,  and  40  A  ,  which  are  equal  to  0.05, 
0.125,  and  0.2,  consecutively,  of  the  far-field  distance  of  each  aperture 
antenna.  The  transmitting  Antenna  A  was  skewed  at  Co  *  -26°  for  all  of  the 
computations  presented  by  these  figures.  Three  curves  are  shown  for  each 


m 

.£%* 

•A 

3S 


$ 

& 


3 

$ 

Oil 


Mm 


■> 

% 


42 


u 

CO  ^ 
H  3 
9  <U 
U  ^ 
M  CO 


a) 

00 

£ 

cO 

u 

a 

o 


tl  flJ  *H 


4J 

CO 


<<  h 

I  o  . 

o<m  ao 

*H  <U 

.*< 

M-t  O 

O  f-t 


■u 

£ 


*4  <*4 

•H  O 
CO 

a  v 
a 

CO  £ 
CO 

U  u 
O  CO 
M-l  *H 
*3 

<D 

tH  e 

oo  o 
c 

CO 


u  a) 
<0  w 
CMC 
O  CO  *H 

*h  a. 
u  <u  ^ 
CO  CO  o 

4J  U-l 

o 
u 


CO  o 
3  4-4 

CO 

^4  CO 

<y  co 

>  c 
u  a) 
a)  a 
3  £ 
O  co 
a. 

0) 

*3  *4 
<u  3 

r-f  U 

a.  V4 

3  <u 

o  a 

CJ  CO 


CO 

0) 

<u 

V4 

00 

a) 

*3 

O  CO 
*  <U 
vo  a 


CSJ  u 

I  00 
<D 

4-i  *3 
O 

o 

3  CN 
rH  rH 
00 

£  «-< 
3  O 


V 

p-4 

00 

ss 


fH 

<44 

3 

01 

0 

Cl 

n 

9) 

01 

•H 

00 

O 

CO 

c 

CO 

n 

u 

1 

o 

0 

<44 

c 

pH 

o 

<44 

,  < 

•H 

m 

44 

O 

CM 

CO 

u 

u 

<44 

00 

•H 

Q 

01 

CO 

44 

a 

ai 

e 

o 

•H 

C0 

e 

to 

C 

n 

44 

CO 

0 

0) 

<44 

•H 

n 

o> 

> 

o 


•n 

0) 

■H  c 
00  o 
C  s 

«  u  O 

CO  «h 

C  U  44 

O  (0  eo 
•nan 
h  <y  oo 
<0  03  4) 

44  44 

o  to  d 

U  *H 

n 

09  o 
3  <44 
03 

n  os 

<U  CO 

>  c 
c 

i-i  ai 
<y  44 

s  c  _ 

O  nj  "O  „ 

a  a) 

VO  V 
-o  u  .  n 
ai  3  vo  oc 

H  U  CN  0) 

a  u  \  *o 

3  a> 

o  a  u*  oo 
O  c0  o  >cr 


u 

o 


03 

ai 

aj 

u 

00 

a) 


09- 


< 


O 

CO 

LU 

111 

a: 

oo 

row 

o 

UJ 

-J 


01 


00 

n 


w 

Ij 

(0 

r-H 

> 

44 

3 

0> 

O 

0 

M 

co 

01 

*H 

00 

O 

C0 

c 

<0 

•< 

M 

M 

1 

O 

o 

44 

a 

rH 

o 

<< 

•H 

4-1 

o 

U 

0 

MT 

c0 

u 

u 

44 

00 

v4 

0 

01 

<0 

11 

CU 

0) 

d 

o 

cO 

c: 

«o 

d 

M 

u 

co 

O 

CO 

44 

t4 

V4 

01 

01 

> 

•H 

c 

0 

o 

c 

c 

to 

_> 

0 

CO 

*H 

c 

M 

u 

o 

CO 

CO 

•w 

3. 

V4 

u 

0) 

00 

CO 

0) 

0) 

u 

14 

o 

c0 

c 

u 

M 

cn 

0 

u 

3 

44 

0 

Cfl 

44 

u 

Cfl 

01 

cO 

CO 

> 

d 

01 

d 

01 

u 

0) 

u 

0J 

■u 

CO 

s 

d 

01 

0 

c0 

T3 

a. 

01 

O 

tj 

u 

• 

01 

d 

vO 

r-4 

•M 

CM 

a 

M 

1 

3 

ai 

O 

O-  44  ■ 

o 

(0 

O  ' 

i 

<r 


29  degrees. 


Figure  4-12.  Coupled  power  versus  rotation  angle  for  a  pair  of  10-X  cir< 
aperture  antennas  for  a  separation  distance  of  40X  for  a  si 
angle  of  -26.0  degrees  for  integration  over  an  integration 


separation  distance.  Each  of  the  three  curves  for  each  separation  distance 
was  computed  for  the  integration  ranges  specified  in  the  caption. 

Figures  4-4,  4-5,  and  4-6  show  the  coupled  power  versus  rotation  angle 
for  the  10-  X  separation  distance  for  integration  over  angular  regions  of 
total  width  70  degrees,  90  degrees  and  120  degress,  consecutively.  The  90 
degree  angle  corresponds  to  the  angle  ^gub  mutually  subtended  by  this  pair 
of  antennas  separated  by  10  X ,  as  computed  from  the  formula 

■®Wb  •  2  arctan 


Pa+Pb 

.  2R  (4-1) 


where  all  symbols  are  as  previously  defined.  The  integration  region  of  70 
degrees  was  chosen  to  investigate  the  impact  on  accuracy  of  using  an 
integration  region  that  extends  beyond  the  mainbeam  but  that  is  smaller 
than  the  mutually  subtended  angle.  The  integration  range  of  120  degrees 
yields  essentially  the  same  coupled  power  curve  as  the  one  obtained  by 
integrating  over  the  entire  front  hemisphere  as  required  by  Equation  (3-7). 
Inspection  of  these  three  curves  shows  that  all  three  curves  have  their 
peak  value  occurring  very  near  Che  rotation  angle  of  g  »  26°  where  the 
apertures  of  Che  two  antennas  are  parallel.  The  computed  peak  values  occur 
at  24.5°,  25.0°,  25.0°  and  for  Che  integration  ranges  of  70  degrees,  90 
deg-.-ees ,  and  120  degrees,  consecutively.  The  curves  for  the  integration 
ranges  of  70  degrees  and  90  degrees  differ  noticeably,  but  not  drastically, 
from  Che  curve  for  the  120  degree  integration  range.  These  limited 
results  indicate  that  the  integration  range  specified  by  Equation  (3-7) 
yields  a  conservative  guideline  provided  that  the  skew  angle  is  much 
smaller  than  the  angular  half-width  of  the  stationary  phase  region.  The 
curves  computed  from  the  integration  ranges  of  70  degrees  or  90  degrees  may 
be  accurate  enough  for  some  EMC  applications. 

The  curves  presented  in  Figures  4-7  through  4-9  for  a  separation 
distance  of  25  X  and  those  presented  in  Figures  4-10  through  4-12  for  a 
separation  distance  of  40 X  further  illustrate  general  trends  and  provide 
additional  understanding  of  the  impact  of  the  integration  region  on 
computational  accuracy.  One  notes  that  the  peak  value  of  coupled  power  has 
shifted  progressively  closer  to  the  stationary  phase  direction,  which 
coincides  with  a  rotation  angle  of  0  ■  0  degrees  on  the  plots,  as  the 


separation  distance  has  been  increased.  In  particular ,  the  peak  value  £or 
the  aost  accurately  computed  curve  (i.e.,  the  one  obtained  with  the  largest 
integration  region)  for  each  separation  distance  occurs  at  an  angle  of 
about  25  degrees,  12  degrees,  and  -3  degrees  for  the  separation  distances 
of  10X  ,  25X  ,  and  40X  ,  consecutively.  This  is  a  sensible  trend  which 
implies  the  well  known  result  that  for  very  large  separation  distances  the 
power  coupled  to  the  receiving  antenna  is  maxisnan  when  its  aainbeam  is 
pointed  directly  at  the  transmitting  antenna. 

A  significantly  sensitivity  in  the  computed  curves  to  the  size  of  the 
integration  range  may  be  discerned  by  inspecting  the  curves  for  the 
separation  distances  of  25X  and  40X  .  The  integration  range  of  48  degreed 
determined  from  Equation  (3-7)  for  the  separation  distance  of  25x  is  also 
approximately  the  angle  mutually  subtended  by  these  two  antennas  for  that 
separation  distance.  Similarly,  the  integration  range  of  29  degrees  for 
the  separation  distance  of  40X  is  also  approximately  the  angle  mutually 
subtended  by  the  two  antennas  for  that  separation  distance.  The 
integration  range  of  70  degrees  provides  integration  over  a  cone  of  half 
angle  equal  to  35  degrees  which  extends  considerably  beyond  the  peak  of  the 
skewed  transmitting  (far-field)  mainbeam.  The  integration  region  of  90 
degrees  yields  coupled  power  plots  that  have  essentially  converged  to  a 
stable  result.  Thus,  it  was  necessary  to  integrate  over  a  range  almost 
twice  the  size  of  the  stationary  phase  region  for  the  separation  distance 
of  25 X  and  to  integrate  over  a  range  almost  three  times  the  size  of  the 
stationary  phase  region  for  the  separation  distance  of  40X  when  the 
transmitting  antenna  is  skewed  at  £0  *  -26  degrees. 

These  results  provide  useful  insight  into  general  coupling  trends  and 
the  impact  of  the  size  of  the  integration  region  on  computational  accuracy. 
However,  it  should  be  recognized  that  some  aspects  of  both  the  general 
trends  and  the  sensitivity  of  the  computational  accuracy  to  the  integration 
range  may  depend  significantly  on  the  aperture  distribution  and,  hence,  on 
the  far-field  pattern.  Accordingly,  it  would  be  useful  to  conduct  similar 
numerical  studies  involving  other  pattern  functions. 

C,  GTSAF  Coupling  Calculations  and  t,.  perimental  Data 

A  sample  execution  for  the  antenna  coupling  algorithm  GTSAF  is 
presented  in  Figure  4-13.  GTSAF  is  designed  to  accept  complex  input 
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"PROGRAM  GTSAF  COMPUTES  THE  NEAR  FIELD  MUTUAL" 

"GAIN  FOR  A  PAIR  OF  ANTENNAS  OPERATING  EITHER  ON" 

"A  CLEAR  SITE  OR  IN  THE  PRESENCE  OF  A  TALL  ROUND  MAST  OR  SHEET" 
"OBSTACLE.  THE  PROGRAMS  EHPARAB  AND  GTOBS  CAN  BE  USED  TO" 
"GENERATE  THE  INPUT  ANTENNA  AND  THE  ANTENNA/ SCATTERER" 

"DATA  FILES  REQUIRED  BY  GTSAF  " 

"LINEAL  DIMENSIONS  ARE  IN  FEET" 

"ENTER  DIAMETER  OF  CIRCUMSCRIBING  SPHERE  FOR  ANTENNA  A" 

? 

4. 

"ENTER  DIAMETER  OF  CIRCUMSCRIBING  SPHERE  FOR  ANTENNA  B" 

? 

4. 

"ENTER  FREQUENCY  IN  GHZ" 

? 

5.5 

"ENTER  THE  DISTANCE  BETWEEN  THE  ANTENNAS" 

? 

16.8 

"ENTER  ANGULAR  COORDINATES  OF  ANTENNA  A  WITH  RESPECT" 

"TO  ANTENNA  B" 

? 

90. 

? 

0. 

"ENTER  THE  AZIMUTH  OFFSET  ANGLE  FOR  ANTENNA  A" 

? 

0. 

"ENTER  ANTENNA  B  AZIMUTH  ROTATION  ANGLE  LIMITS" 

? 

-24. 

? 

24. 

"ENTER  ANGULAR  INCREMENTS  FOR  ROTATION  ANGLE" 

7 

1. 

"END  OF  EXECUTION.  OUTPUT  MUTUAL  GAIN  DATA"’ 

"IS  IN  FILE  MGOUT . " 
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Figure  4-13.  Sample  execution  for  the  antenna  coupling  computer 
program  GTSAF  for  the  computed  data  plotted  in 
Figure  4-15. 
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ROTATION  ANGLE 

Figure  4-14  .  Antenna  coupling  versus  antenna  rotation  angle  for 
the  near-field  situation  depicted  in  the  inset  of 
Figure  4- 15  with  the  round  mast  obstacle  removed, 
i.e.,  on  a  "clear  site." 
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ROTATION  ANGLE 

Figure  4-15  .  Antenna  coupling  versus  antenna  rotation  angle 
for  the  near-field  situation  depicted  in  the 
inset  for  5.5  GHz. 


antenna  data  and/or  composite  antenna/obstacle  data  that  has  been  measured 
or  data  that  has  been  generated  from  an  auxiliary  program.  The  GTSAF 
algorithm  can  be  employed  to  compute  near-field  power  coupling  for  pairs  of 
reflector  antennas ,  wire  antennas ,  horn  antennas ,  and  phased  arrays . 
Coupling  between  these  different  types  of  antennas  can  be  based  on  a 
knowledge  of  their  SAF's.  The  GTSAF  results  for  two  reflector  antennas 
operating  on  a  clear-site  and  in  the  presence  of  a  round  mast  obstacle  are 
shown  in  Figures  4-14  and  4-15,  respectively.  The  measured  coupling  data 
.or  Figures  4-14  and  4-15  were  obtained  during  a  previous  Navy  contract. 

The  GTSAF  coupling  data  for  these  two  figures  were  obtained  by  using 
theoretical  SAF's,  as  follows.  First,  the  obstacle  effects  program  GTOBS 
was  employed  to  generate  both  clear-site  and  blocked  pattern  data.  Cosine- 
on-alO-dB  pedestal  aperture  illumination  functions  were  used  for  the  two 
reflectors.  The  GTOBS  computer  SAF  data  were  then  used  as  input  data  to 
GTSAF  to  compute  the  coupling  data  displayed  in  Figures  4-14  and  4-15. 
Inspection  of  the  antenna  coupling  for  the  blocked  situation  shown  in 
Figure  4-15  reveals  that  a  good  engineering  description  of  the  measured 
coupled  power  is  obtained  from  GTSAF  despite  the  fact  that  the  computed 
unobstructed  coupled  power  data  do  not  closely  replicate  the  "sidelobes"  of 
the  measured  unobstructed  power  coupling  data  displayed  'n  Figure  4-14. 

D,  Auxiliary  Programs  GTOBS  and  EHPARAB 

Program  GTOBS  is  designed  to  accept  antenna  patterns  generated  from 
EHPARAB  and  to  complete  blocked  patterns  using  either  cylinders  or  flat 
plates  as  the  scattering  obstacles  located  near  the  antenna.  These 
patterns  then  can  be  written  to  files  for  use  as  input  to  the  mutual 
coupling  program  GTSAF.  Initial  input  to  GTOBS  is  partially  from  files  and 
partially  interactive.  It  is  assumed  that  before  program  execution 
elevation  and  azimuthal  antenna  patterns  are  contained  in  complex  form  in 
files  ANTEL  ard  ANTAZ,  respectively.  Each  complex  value  is  preceded  by  the 
corresponding  theta  or  phi  angle.  The  first  column  of  the  elevation  file 
contains  an  elevation  cut  for  the  front  hemisphere,  i.e.,  <t>  *  0°, 

0°<  9<180°.  The  second  column  contains  a  sir  -i ir  cut  for  the  back 
hemisphere,  for  <p  m  180°,  0°<9<180o.  The  azimuthal  file  contains  an 
azimuth  cut  for  the  entire  9  »  90°  plane,  0°<  $<360°.  The  initial 
interactive  input  is  mostly  concerned  with  the  physical  parameters  of  the 


situation  for  which  a  pattern  ia  to  be  computed.  Information  auch  aa  the 
width  of  the  antenna,  the  position  of  the  obatacle  relative  to  the  antenna, 
and  the  shape  and  size  of  the  obstacle  is  input  in  this  section. 

EHPARAB  is  used  to  compute  the  E-plane  and  H-plane  patterns  of  a  horn- 
fed  parabolic  dish  antenna  over  complete  circles  encompasing  the  antenna. 
EHPARAB  is  based  on  the  "physical  optica"  integral  over  the  reflector 
current  distribution,  as  described  by  Silver  in  Microwave  Antenna  Theory 
and  Design.  The  computed  complex  electric  field  data  are  then  used  in 
GTOBS  to  compute  composite  antenna/scat terer  patterns  and  in  GTSAF  to 
compute  antenna  coupling.  Plots  of  the  computed  amplitude  of  the  electric 
field  scattered  from  the  reflector  dish  due  to  illumination  by  a  short 
dipole  located  at  the  focus  are  shown  in  Figure  4-16  and  4-17  for  the  E-and 
H-planes,  respectively.  The  EHPARAB  data  precisely  match  other  published 
data  for  this  case.  The  operator  inputs  to  EHPARAB  are  frequency,  the  F/D 
ratio  for  the  dish,  the  peak  gain  and  the  principal  plane  pattern  functions 
of  the  feed  horn,  and  the  direction  of  the  assumed  linear  polarization  of 
the  feed  horn. 
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Figure  4-16.  E-plane  scattered  pattern  for  a  reflector 
antenna  fed  by  a  short  dipole  for  F/D  =  0 
D  =  10  wavelengths. 
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SECTION  V 

CONCLUDING  KSMABKS  AND  RECOMMENDATIONS 


The  efforts  to  develop  an  efficient  and  versatile  computer  algorithm 
for  computing  antenna  coupling  in  the  presence  of  near-field  obstacles  were 
successful.  The  results  of  the  analytical  and  numerical  research  work  show 
that  the  algorithm  based  on  the  Spherical  Angular  Function  (SAF)  technique 
can  be  used  for  applications  involving  antennas  of  low,  moderate  and/or 
high  directivity. 

The  basic  SAF  equations  were  first  derived  for  describing  clear-site 
antenna  coupling  in  terms  of  the  normalized  SAP's  of  the  antennas.  The 
normalized  SAF  F  ( 0 , )  of  an  antenna  is  related  to  the  far-field  complex 
electric  field  E^Ce,^)  of  the  antenna  as  Fq(6,$)  *  [g]  reJ  Ft^(0,4>), 
where  r  is  the  far-field  range  of  the  observation  point  and  g  is  the  peak 
far-field  gain  of  the  antenna.  The  coupling  analysis  requires  the 
evaluation  of  double  convolution  integrals  involving  the  SAF's  of  the  two 
antennas  over  portions  of  the  far-field  sphere,  where  the  notion  of 
convolution  was  extended  to  include  the  type  of  integrals  encountered  in 
near-field  antenna  coupling  analysis.  The  corresponding  SAF  equations  for 
antennas  operating  in  the  presence  of  near-field  scattering  obstacles  were 
also  derived. 

Numerical  data  are  calculated  and  studied  to  determine  the  general 
capabilities  of  the  SAF-based  antenna  coupling  algorithm  and  two  related 
auxiliary  algorithms.  Firstly,  coupling  data  were  computed  for  a  pair  of 
*S-A  dipoles  and  for  a  pair  of  10-A  circular  aperture  antennas.  Secondly, 
mutual  coupling  data  were  calculated  for  a  pair  of  directive  antennas  for 
clear-site  as  well  as  blocked  conditions  in  the  near-field.  Good  agreement 
was  obtained  between  measured  and  computed  data.  The  coupling  program  is 
denoted  as  GTSAF,  and  is  supplied  to  CECOM  under  separate  transmittal.  The 
auxiliary  program  GTOBS  and  EHPARAB,  for  computing  composite 
antenna/obstacle  SAF's  and  for  generating  SAF's  for  paraboloidal  reflector 
antennas,  respectively,  were  also  implemented  and  validated. 

It  would  be  particularly  beneficial  to  CECOM  for  Georgia  Tech  (1)  to 
make  certain  refinemets  and  extensions  of  the  coupling  prediction 
capabilities  developed  under  this  contract  and  (2)  to  obtain  additional 


measured  dace  for  model  validations.  Accordingly,  Che  following  research 
tasks  are  proposed  in  order  to  provide  additional,  validated  analytical 
capabilities  to  CECOM  for  monitoring  and  improving  EMC: 

Task  1.  Rapid  Determination  of  SAP* a  for  Antennas,  Scatterera,  and 
Equipments .  The  SAF's  of  antennas,  scatterers,  or  equipments  may  be 
readily  computed  from  the  complex  field  data  obtained  from  Near-field  (NF) 
measurements.  The  feasibility  of  exploiting  the  Modulated  Scattering 
Technique  (MST)  to  develop  a  portable  near-field  test  facility,  capable  of 
acquiring  the  NF  data  in  a  few  minutes  or  seconds,  should  be  investigated. 
A  rapid  NF  measurement  capability  combined  with  the  GTSAF  algorithm  will 
significantly  enhance  CECOM' a  ability  to  determine  and  control  EMI. 

Task  2.  Out-of-Band  Computer  Model  Simulations.  A  Monte  Carlo 
computer  model  has  been  developed  under  a  previous  government  contract  to 
compute  statistical  average  out-of-band  reflector  antenna  patterns. 
Preliminary  computer  simulations  have  demonstrated  that  the  out-of-band 
patterns  are  particularly  sensitive  to  the  relative  phase  deviation  of  the 
higher-order  transmission  line  mode.  However,  the  sensitivity  of  the  model 
results  as  a  function  of  variables  such  as  the  mode  amplitudes  mode  phases, 
and  "mixture"  of  modes  has  not  yet  been  tested  against  measured  data.  It 
is  proposed  that  the  following  numerical  simulations  be  performed  to 
determine  the  computer  model  sensitivity: 

a.  Exercise  the  computer  model  to  determine  the  effects  of  mode 
amplitudes,  mode  relativae  phases,  and  mode  mixture  on  the  out- 
of-band  pattern  statistics  --  specifically,  the  mean  pattern  and 
standard  deviations. 

b.  Modify  the  model  to  account  for  the  effects  of  waveguide  feed 
system  length  on  the  relative  model  phases  by  using  the  higher- 
order  mode  phase  velocity  equations. 

c.  Compare  the  computer  model  mean  patterns  and  standard  deviations 
with  the  measured  data  obtained  during  the  measurement  program. 

Task  3.  Extension  of  the  Out-of-Band  Reflector  Antenna  Pattern  Model. 
The  present  computer  model  performs  a  vector  calculation  of  the  average 


64 


pattern  and  standard  deviation  over  an  approximately  +50°  to  ±80°  field-of- 
view  centered  on  the  antenna  boresight.  In  many  applications,  it  is 

desirable  to  determine  the  complete  pattern  statistics  including  the  back- 
lobe  regions.  Also,  it  is  often  necessary  to  compute  pattern  statistics 
such  as  the  median  gain  and  standard  deviation  over  either  a  complete  4tt 
steradian  or  partial  angular  sector.  However,  with  the  present  model, 
these  computations  can  be  time-consuming  and  expensive.  Thus,  it  is 
proposed  that  the  model  be  extended  to: 

a.  Include  the  effects  of  reflector  rim  edge  diffraction  and  strut 
and  feed  horn  blockage  on  the  out-of-band  antenna  pattern.  This 
extension  will  make  it  possible  to  compute  the  out-of-band 
pattern  over  a  full  4ir  steradian  angular  sector,  or  over  any 
selected  subsector  of  the  pattern. 

b.  Incorporate  routines  to  determine  the  median  gain  and  standard 
deviation  of  the  out-of-band  pattern  for  the  full  4ir  steradian 
angular  sector. 

c.  Examine  computational  methods  to  decrease  the  computer  run  time 
required  by  the  model.  These  methods  might  incorporate  some 
recently  published  techniques  which  are  particularly  suited  to 
circular  aperture-type  antennas. 

Task  4.  Additional  Out-of-Band  Pattern  Measurements.  Additional 
measurements  of  the  out-of-band  pattern  performance  of  reflector  antennas 
are  needed  to  determine  the  effects  of  transmission  line  components.  The 
out-of-band  pattern  performance  has  been  shown  to  be  sensitive  to  the  type 
of  waveguide  transmission  components  employed  in  the  feed  system.  To  date, 
measurements  have  been  performed  using  a  simple  coaxial  cable-to-waveguide 
adapter  feed  system.  It  is  proposed  that  additional  measurements  be 
performed  using  a  feed  system  incorporating: 

a.  coaxial  cable-to-waveguide  adapter  and  E-plane  bend, 

b.  coaxial  cable-to-waveguide  adapter  and  H-plane  bend, 

c.  a  coaxial  cable-to-waveguide  adapter  attached  to  a  typical 
waveguide  feed  which  extends  from  the  reflector  rim  to  the  feed 
horn. 
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Statistical  processing  will  be  performed  on  the  data  to  determine 
whether  the  different  transmission  line  devices  result  in  significant 
changes  in  the  out-of-band  statistical  average  pattern  performance.  These 
measurements  will  be  performed  over  frequency  bands  centered  about  the 
fundamental,  the  second,  and  the  third  harmonic  frequencies. 

Task  5.  Out-of-Band  Antenna  Coupling.  Very  little  measured  data  are 
available  validate  out-of-band  data  computed  via  GTSAF.  It  is,  therefore, 
proposed  to  perform  a  limited  set  of  out-of-band  antenna  coupling 
measurements ,  including  the  following: 

a.  Out-of-band  coupling  measurements  as  a  function  of  azimuth  angle 
for  two  reflector  anten.as  located  at  the  same  elevation.  (Note: 
this  corresponds  to  the  “common  boresight"  case.)  These 
measurements  will  be  p;-  rformed  over  frequency  bands  centered  on 
the  second  and  third  ha?  .Tonics . 

b.  Antenna  coupling  measure  ts  as  a  function  of  azimuth  angle  for 
two  reflector  antennas  separated  by  approximately  one-to-two 
antenna  diameters  in  elevation.  These  measurements  will  be 
performed  over  frequency  bands  centered  on  the  second  and  third 
harmonics . 

The  measured  data  will  be  processed  to  determine  the  mean  and  standard 
deviation  of  the  out-of-band  coupling  as  a  function  of  azimuth  angle. 


snrnoi  vi 


1KRBEHCBS 


1.  R.  F.  Harrington,  Field  Computation  by  Moment  Methods,  MacMillan  Co., 

Nev  York  1968.  ~ 

2.  A.  T.  Adams,  "Method  of  Moments  Applications-Vol .  I.,  An  Introduction 
to  the  Method  of  Moments,"  P.ADC-TR-73-217,  November  1974,  AD  #A- 
002820. 

3.  R.  C.  Hansen,  Microwave  Scanning  Antennas,  Volume  I,  Academic  Press, 
New  York,  1964. 

4.  F.  L.  Cain,  C.  E.  Ryan,  B.  J.  Cown,  and  E.  E.  Weaver,  "Electromagnetic 
Effectiveness  Investigations  of  Near-Field  Obstacle  Effects,  Antenna 
Coupling  and  Phased  Arrays,"  Georgia  Institute  of  Technology,  Final 
Engineering  Report,  Contract  No.  N00024-72-C-1276. 

5.  R.  C.  Rudduck  and  R.  J.  Markefka,  "Quasi-Optical  Techniques  for 
Antennas  at  UHF  and  Above,"  Ohio  State  University  ElectroScience 
Laboratory,  Quarterly  Report  712242-6,  Contract  No.  N00123-79-C-1469, 
February  1981. 

6.  B.  Ye.  Rinber,  "Decoupling  of  Closely  Situated  Reflector  Aerials," 
Radiotekhnika  i  Electronika  6,  No.  6,  1961,  pp.  907-916. 

7.  E.  Schweicher,  "Linkage  Between  Two  Microwave  Antennas  in  the  Presence 
of  a  Cylindrical  Reflector,"  Annales  Des  Telecommunications,  Tome  34, 
Nos.  9-10,  September-October  1979. 

8.  C.  E.  Ryan,  Jr.,  E.  E.  Weaver,  and  W.  P.  Cooke,  "Development  of  Total 
Electromagnetic  Effectiveness  During  Ship  Design,"  Georgia  Institute 
of  Technology,  Final  Technical  Report,  Contract  No.  N00024-77-C-5338, 
December  1979. 

9.  C.  E.  Ryan,  Jr.,  E.  E.  Weaver,  and  B.  J.  Cown,  "Plane  Wave  Spectrum 
Analysis  of  Near-Field  Obstacle  Effects  on  Directive  Antenna 
Patterns,"  IEEE  Transactions  on  Antennas  and  Propagation,  27,  No.  6, 
November  1979. 

10.  B.  J.  Cown  and  C.  E.  Ryan,  Jr.,  "Near-Field  Theory  and  Techniques 

for  Wideband  Radiating  System  at  In-Band  and  Out-of-Band  Frequencies," 
Georgia  Institute  of  Technology,  Final  Technical  Report,  Contract 
No.  DAAG29-78-C-0029,  January  1982. 

11.  A.  D.  Yaghjian,  "Efficient  Computation  of  Antenna  Coupling  and  Fields 
Within  the  Near-Field  Region,"  IEEE  Transactions  on  Antennas  and 
Propagation,  Volume  AP-30,  No.  1,  pps.  113-128,  January  1982. 


67 


i 


1 

*• 

S 

«s 


12.  B.  J.  Corn,  "In-Band  Antenna  Coupling  Research,"  Georgia  Institute 
o£  Technology,  Interim  Technical  Report  No.  1,  Contract  DAAB07-82- 
C-J069,  September  30,  1982. 

13.  G.  D.  Ber gland,  "A  Guided  Tour  of  the  Fast  Fourier  Transform," 

IEEE  Spectrum,  July  1969. 

14.  D.  M.  Kerns,  Plane  Wave  Scattering  Matrix  Theory  of  Antenna  and 
Antenna-Antenna  Interactions,  0.  S.  Government  Printing  Office, 
Library  of  Congress  Catalog  Card  Numbers  81-600015,  June  1981. 

15.  V.  H.  Rumaey,  "Reaction  Concept  in  Electromagnetic  Theory,"  Physical 
Review,  Volume  94,  No.  6,  pp.  1483-1491,  June  15,  1954. 


1  16.  R.  F.  Harrington  and  A.  T.  Villeneuve,  "Reciprocity  Relationships 

1  for  Gyrotropic  Media,"  IRE  Transactions  on  Microwave  Theory  and 

*  Techniques .  MTT-6,  pps.  308-310,  July  1958. 

» 

J  17.  Microwave  Antenna  Theory  and  Design,  S.  Silver  (Editor),  Dover 

*  Publications,  New  York,  1965. 

I 

'l 

)18.  P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics,  Volume  2, 

McGraw-Hill  Book  Co.,  New  York,  1953. 


«  19.  A.  Banos,  Jr.,  Dipole  Radiation  in  the  Presence  of  a  Conducting  Half- 

J  Space,  Pergamon  Press,  New  York,  1966. 

I  ' 

I  20.  L.  Schwartz,  "Theorie  des  Distributions  a  Valeura  Vectorielles ,” 

J  Anmales  de  L'lnstitut  Fourier,  Universite  de  Grenoble,  Tome  8-9, 

J-  ,  1958-1959,  pp.  1-202. 

s 

*  21.  L.  Schwartz,  Mathematics  for  the  Physical  Sciences,  Addison-Wesley 

\  Publishing  Co.,  Reading,  Massachusetts,  1966. 

,  22.  N.  Zittron  and  S.  N.  Karp,  "Higher  Order  Approximations  in  Multiple 

\  Scattering  II.  Three  Dimensional  Scalar  Case,"  Journal  of  Mathe- 

•  matical  Physics,  Vol.  2,  No.  3,  May- June  1961. 

» 

i  23.  K.  Hongo,  "Multiple  Scattering  by  Two  Conducting  Circular  Cylinders," 

j  IEEE  Antennas  and  Propagation,  Vol.  AP-26,  No.  5,  September  1978. 


j  24.  P.  Debeye,  "Approximate  Formulae  for  the  Cylinder  Functions  for 

<  Large  Values  of  the  Argument  and  Any  Value  of  the  Index,"  The 

1  Collected  Papers  of  Peter  J.  W.  Debeye,  Interscience  Publishers, 

I  Inc.,  New  York,  1954. 

'  25.  P.  M.  Morse  and  H.  Feschback,  Methods  of  Theoretical  Physics, 

.  Volume  1,  McGraw-Hill  Book  Co.,  New  York,  1953. 


I 

/ 

i  . 


/  / 


.*  /  v  . 
\  *,  *, 

*  *  »  m 


68 


APPENDIX  A 


DISTRIBUTION  LIST 


i 


/ 


DISTRIBUTION  LIST 


Address  No.  of  Copies 

*  Defense  Technical  Information  Center  * 

ATTN:  DT1C-FDAC 
Cameron  Station  (81dg  5) 

Alexandria,  VA  22304-6145 

Director  01 

National  Security  Agency 
ATTN:  TDL 

Fort  George  G.  Meade,  MO  20755 

DCA  Defense  Communi cat ions  Engineering  Center  01 

COOE  R123,  Tech  Library 
1860  Wiehle  Avenue 
Res  ton,  VA  20090 

Defense  Communications  Agency  •  02 

Technical  Library  Center 
CODE  205  (P.A.  TOLOVI) 

Washington,  OC  20305 

Commanding  Officer  01 

Naval  Research  Laboratory 
ATTN:  COOE  2627 
Washington,  DC  20375 

Commander  01 

Naval  Ocean  Systems  Center 

ATTN:  Library 

San  Diego,  CA  92152 

Commander  01 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Library,  CODE  WX-21 

Commandant  02 

Marine  Coros 
ATTN:  CODE  LMC 
Washington,  DC  20380 

Headquarters  01 

US  Marine  Corps 
ATTN:  CODE  I  NTS 
Washington,  OC  20380 

(♦Note  to  contractor:  These  2  copies  will  be  sent  from  the  STINFO  Office, 

Fort  Monmouth,  N.J. ) 


.'■A'.'.'jV/- 


D I ST RI8UTI0N  LIST  (CONTO) 


Address 

Co, irr, and.  Control  4  Communications  Division 
Development  Center  . 

Marine  Corps  Development  4  Education  Command 
Quantico,  VA  22134 

Naval  Teleccimiunications  Command 
Technical  Library,  CODE  9IL 
4401  Massachusetts  Avenue,  NW 
Washington,  DC  20390 

Naval  Air  Systems  Command 
COOE:  AIR-5332 
Washington,  OC  20360 

Dr.  T.  G.  Berl incourt 

Office  of  Naval  Research  (CODE  420) 

800  N.  Quincy  Street 
Arl ington,  VA  22217 

AUl/LSE  64-285 
Maxwell  AFB,  AL  36112 

Ron*  Air  Development  Center 
ATTN:  Documents  Library  (TSLO) 

Griffiss  ARB,  NY  13441 

USAFETAC/C8TL 
ATTN:  Librarian 
STOP  825 

Scot  AFB,  IL  62225 

Air  Force  Geophysics  Lab 
L.  G.  Hanscom  AFB 
ATTN:  LIR 
Bedford,  MA  01730 

AFGL/SULL 

S-29 

HAFB,  MA  01731 

Keadouarte rs 
AFEWC 
ATTN:  EST 

San  Antonio,  TX  78243 


No,  of  Cooies 
01 


01 


04 


01 


01 

01 


01 


01 


01 


02 


j £-t 


/ 


■>  * 
$ 

m 


m*  if 


DISTRIBUTION  LIST  (CONTD) 


Address 


No.  of  Conies 


Headquarters 

Air  Force  Systems  Command 
ATTN:  DLWA/Mr.  P.  Sandler 
Andres  AFB 

Washington,  DC  20331 

Commander 
MICOM . 

Redstone  Scientific  Information  Center 
ATTN:  Chief,  Document  Section 
Redstone  Arsenal ,  AL  35309 

Commander 

MICOM 

ATTN:  DRSMI-RE  (Mr.  Pittman) 

Redstone  Arsenal ,  AL  35809 

Commandant 

US  Army  Aviation  Center 
ATTN:  ATZq-D-MA 
Fort  Rucker,  AL  36362 

Commander 

Headquarters 

Port  Huachuca  ,  „  .  . 

ATTN:  Technical  Rererence  Division 

Fort  Huachuca,  AZ  85613 
Commander 

US  Army  Electronic  Proving  Ground 

ATTN:  STEEP-MT 

Fort  Huachuca,  AZ  85613 

Commander 

US  Army  Proving  Ground 

ATTN:  STEYP-MTD  (Tech  Library) 

Yuma,  AZ  85364 

Di  rector 

US  Army  Air  Mobility  R&D  Lab 
ATTN:  T.  Gossett,  Bldg  207-5 
Nasa  Ames  Research  Center 
Moffett  Field,  CA  94035 

HQOA  (OAMO-TCE) 

Washington,  DC  20310 


01 


02 


01 


01 


01 


01 


01 


01 


02 


& 


•■.V, 


DISTRIBUTION  LIST  (CONTD) 


No.  of  Copies 


Address 


Deputy  For  Science  &  Technology 
Office,  Asst  Sec  Army  (R&D) 

Washington,  DC  20310 

HQDA  (DAMA-ARZ-O/Or.  F.  D.  Verderame) 
Washington,  DC  20310 

Commandant 

US  Army  Signal  School 
ATTN:  ATSH-CD-MS-E 
Fort  Gordon,  GA  30905 

Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-MS-E 
Fort  Senning,  GA  31905 

Commander 

ARPXOM  •  . 

ATTN:  Systems  Analysis  Office,  DRSAR-PE 
Rock  Island,  IL  61299 

Director,  Combat  Developments 
US  Army  Armor  Center 
ATTN:  ATZK-CO-MS 
Fort  Knox,  KY  40121 

Commander 

US  Army  Test  &  Evaluation  Command 

ATTN:  DRSTE-CT-C 

Aberdeen  Proving  Ground,  MO  21005 

Commander 

Harry  Diamond  Laboratories 
ATTN:  Library 
2800  Powder  Mill  Road 
Adel  phi,  MO  20783 

Oi  rector 

US  Army  Ballistic  Research  Labs 
ATTN:  DRX3R-L8 

Aberdeen  Proving  Ground,  MO  21005 


02 


01 

01 


01 


01 


02 


.  01 


01 


01 


n? ATRlBUTIQN  LIST  (CONTDl 


Address 


Director  .  .  ,  . 

JS  Army  Human  Engir iesr’n9  ^ 
Aberdeen  Proving  Ground,  MO  Zluusj 


US  Material  Systems  Analysis  Activity 
Aberdeen  Proving  Ground,  MO  21005 


SlTr^tiaterlal  Sy»t«M  Analysis  Activity 
ATTN:  ORXSY-MP 

Aberdeen  Proving  Ground,  MO  210US 

CECOM  Aviation  Electronics  Office 
ATTN:  DRSEL-MME-LAF(Z) 

St.  Louis,  MO  63166 


Commander 

AVSCOM 

ATTN:  ORSAV-E 
PO  Sox  209 
St.  Louis,  MO 


Commander 

ATTN:  DRDAR-L1N-S  (Bldg  95} 
Oover,  MJ  07801 


Colander 

ARRAOCOM 

ATTN:  OROAf-TSS,  #59 
Oover,  NJ  07801 


Director  .  ,Ta, 

Joint  Co~unication  Office  L « 
ATTN:  TT-AO  (Tech  Ooc  Cen). 
Fort  Monmouth,  NJ  07703 


TAC) 


PM,  F 1 REF  I ND£ R/ REM3ASS 

ATTN:  DRCPM-FER 

Fort  Monmouth,  NJ  07703 


Corwiande  r 
CECOM 


ATTN: 
ATTN: 
Ft.  Mo 


OELSD-L 

oelsD-l-s(stinfo) 

nnouw,  NJ  07703  '5301 


f 


No.  of  Copies 
01 

01 


01 


01 


01 


01 


01 


01 


01 


This  Document 
Reproduced  From 
“**•  Available  Copy 


DISTRIBUTION  list  (COHTO). 


Address 

Project  Manager 
NAVCON 

ATTN:  ORCPM-NC-TM 
Fort  Monmouth,  NJ 


Commander 

US  Army  Satellite 
ATTN:  ORCPM-SC-3 
Fort  Monmouth,  NJ 


Communications  Agency 
07703 


TRI-TAC  Office 

ATTN:  TT-SE 

Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Avionics  Lab 
AVRADA 

ATTN:  SAYAA-D 

Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Research  Office 
ATTN:  Or.  Horst  Hittmann 
PO  Box  12211 

Research  Triangle  Park,  NO 


27709 


Cormander 

US  Army  Research  Office 
ATTN:  ORXRO-IP 
PO  Box  12211 
Research  Triangle  Park, 


NC 


27709 


Commander  _  . 

US  Army  Research  Office 
ATTN  ORXRC-PK  (Or.  R.  J 
Research  Triangle  Park, 


.  Lontz) 
NC  27709 


Commandant  .  ..  .•_  _ 

US  Army  Inst  For  Military  ^ssis^cnce 

ATTN:  ATSU-CTQ-MO ^ 

Fort  Bragg,  NC  2S307 


Commandant 

U2AFAS 

ATTN:  ATSr-CO-OC 
Fort  Sill,  OK  7 ao03 


No. ‘of  Copies, 
01 


02 


01 


01 


01 


02 


02 


02 


01 


i 


This  Document 
Reproduced  From 
Best  Available  Copy 


DISTRIBUTION  LIST  (CCNTQ) 


Address 

Cc~andant 

US  Army  Air  Defense  School 
ATTN :  ATS A -CD -MS -C 
Fort  81iss,  TX  79S16 


Commander 
AMC  ■ 

ATTN:  AMCDE 

5001  Eisenhower  Avenue 

Alexandria,  VA  22333.0001 

Director 

US  Army  Signals  Warfare  ctr. 

Vint  hill  Farms  Station  ATTN:  AMSEL-SW-OS 
Warrentcn,  VA  22186  -5100 


I  I  CL  VU  I  .  -  ,  .  - 

US  Army  Signals  '.Jarre re  ctr. 
ATTN:  ATTN:  AMSEL-SW-A0 
Vine  Hill  Farms  Station 
Varrenton,  VA  22i8A -5100 


Commandant 

US  Army  Engineer  School 

ATTN:  ATZA-T0L 

Fort  3el voi r,  VA  22C60 

Commander 

US  Army  Engineer  Topographic  Labs 

ATTN:  ETL-TD-EA 

Fort  3elvoir,  VA  22050 

Commander 

US  Army  Logistics  Center 
ATTN:  ATCL-MC 
Fort  Lee,  VA  22391 

Commander 
TP/1  DOC 

ATTN:  TDOC-TA 

Fort  Monroe,  VA  23551 


US  Army  Training  l  Doctrine  Co...  and 

nl  i  N :  AT  i_C -  i  in 

Fort  Monroe ,  VA  2 7S5 1 


Mo.  of  Copies 
01 


01 


02 


01 


02 


01 


02 


01 


01 


This  Document 
Reproduced  From 
Best  Available  Copy 


distribution  list  (conto) 


Address 


fJo. 


Commander 

US  Army  Garrison 

Vint  Hill  Farms  Station 

ATTN:  IAVAAF 

War  ran  ton,  VA  221S6 

Project  Manager 
Control  &  Analysis  Centers 
Vint  Hill  Farms  Station 
Warrenton,  VA  22186 

Commander 

Niqht  Vision  &  Electro-Optics  Lab 
CECOM 

ATTN:  CELMV-0 

Fort  Bel  voir,  VA  22C60 


Commander 

Atmospheric  Sciences  Lab 
LABCOM 

ATTN:  SLCAS-ST_S 

White  Sands  Missile  Range,  NM  83002 
Ch  i  e  f 

Office  Missile  Electronic  Warfare 
Electronic  Warfare  Lab,  LABCOM 
V/hite  Sands  Missile  Range,  NM  83002 

Chief 

Intel  Material  Development  4  Support  Office 
Electronic  Warfare  Lab,  ESAOCOM 
Fort  George  G.  Meade,  MO  20755 

C  C  1  1  i*  2  Pi  Q  F 

US  Army  Ballistic  Research  Lab/ARRADCCM 
ATTN:.  CRCAR-T52-S  (ST  INFO) 

Aberdeen  Proving  Ground,  MD  '  21005 

Co~ande  r 
LABCOM 

ATTN:  AMSLC-CT 
2800  Pc.; r  Mi  1 1  Road 
Ade'p.ni,  '0  207'-?-l’A5 

Co  .-.render 
LABCOM 


of  Conies 
01 


01 


01 


01 


01 


01 


01 


02.  • 


01 


This  Document 
Reproduced  From 
Best  Available  Copy 


DISTRIBUTION  LIST  (CQMTD) 


Address 

Headquarters 

Harry  Piarnond  Laboratories 
ATTN:  SLCHD  -TO  (Or.  U.  W.  Carter) 

2SOO  Powder  Mill  Road  • 

Adel  phi,  MO  20783  -1145 

MIT  -  Lincoln  Laboratory 
ATTN:  Library  (RM  A-082) 

PO  Box  73 

Lexington,  MA  02173 

NASA  Scientific  &  Tech  Information  Facility 
Baltirrore/Washington  Inti  Airport 

PO  Box  8757,  MO  21240 
• 

National  Bureau  of  Standards 
SI dg  225,  RM  A- 331 
ATTN:  Mr.  Leedy 
Washington ,  OC  20231 

Dr.  Clayton  Paul 
Electrical  Engineering  Oept 
University  of  Kentucky 
Lexington,  KY  40506 

Mr.  John  Spina 
RAQC/RBCT 

Griffiss  AFB ,  NY  13441 
Commander 

US  Army  Communications.  Research  A  Oevelooment  Command 
ATTN:  ORCCO-SEI -A/Mr .  Sam  Segner 
Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Cofr.muni  cat  ions -Electronics  Command 
ATTN:  AMSEL  -COM-D 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Communications  -  Electronics  Command 
ATTN:  AMSEL-C0M-AY-2/Mr .  Stuart  AlLert 

Fort  Monmouth ,  NJ  07703-5202 


No.  of  Copies 

x  01 


01 


01 


01 


01 


01 


01 


15 


This  Document 
Reproduced  From 
Best  Available  Copy 


